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ABSTRACT 

We present four new epochs of K s -band images of the young pre-transitional disk around LkCa 15 and perform extensive forward 
modeling to derive the physical parameters of the disk. We find indications of strongly anisotropic scattering (Henyey-Greenstein 
g = 0.67/jJ |j) and a significantly tapered gap edge (“round wall”) but see no evidence that the inner disk, whose existence is predicted 
by the spectral energy distribution, shadows the outer regions of the disk visible in our images. We marginally confirm the existence of 
an offset between the disk center and the star along the line of nodes; however, the magnitude of this offset ( x = 27 /kj mas) is notably 
lower than that found in our earlier H - band images. Intriguingly, we also find an offset of y = 69/ 49 mas perpendicular to the line of 
nodes at high significance. If confirmed by future observations, this would imply a highly elliptical - or otherwise asymmetric - disk 
gap with an effective eccentricity of e ~ 0.3. Such asymmetry would most likely be the result of dynamical sculpting by one or more 
unseen planets in the system. Finally, we find that the bright arc of scattered light we see in direct imaging observations originates 
from the near side of the disk and appears brighter than the far side because of strong forward scattering. 

Key words, stars: individual: LkCa 15 - stars: pre-main sequence - planetary systems - planets and satellites: formation - 
protoplanetary disks 


1. Introduction 

As potential indicators of planetary companions, transitional 
disks offer the tantalizing possibility of observing planet for- 
mation in action. They were initially identified as protoplane- 
tary disks with a reduced near-infrared (NIR) excess, indicating 
depleted inner regions (Strom et al. 1989; Calvet et al. 2005). 
Millimeter-wave imaging has confirmed that these disks indeed 
contain large inner holes or annular gaps (e.g., Andrews et al. 
201 1). However, optical and NIR images do not always reveal 
these holes (Dong et al. 2012), which indicates different spatial 
distributions of micrometer and millimeter-sized dust grains (see 
de Juan Ovelar et al. 2013). In addition, stars continue to accrete 


* Appendices are available in electronic form at 
http : //www. aanda. org 


substantial amounts of gas despite being cut off from the main 
gas reservoir in the outer disk (Ingleby et al. 2013; Bergin et al. 
2004), which increases the complexity of the puzzle. 

To understand the complex geometry of these disks, multi- 
wavelength imaging is necessary. Only a handful of these gaps 
have been imaged in the optical/NIR (e.g., Mayama et al. 2012; 
Hashimoto et al. 2012; Quanz et al. 2013; Avenhaus et al. 
2014; Boccaletti et al. 2013; Garufi et al. 2013), due to their 
high contrast ratio and proximity to the host star. A promising 
method to overcome these hurdles is angular differential imag- 
ing (ADI; Marois et al. 2006), which utilizes the natural field 
rotation of altitude-azimuth ground based telescopes to improve 
high-contrast imaging sensitivity. The ADI method efficiently 
reduces the impact of the stellar point spread function (PSF) 
wings, though at the cost of also imposing a non-negligible 
amount of self-subtraction for off-axis sources, such as planets 
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Table 1. Summary of high-contrast observations of LkCa 15. 


ID 

Epoch 

Band 

DIT (s) 

Co-adds 

# Frames 

Time (min) 

FoV 

Parall. angles (°) 

Sat. 

PSF stars 

Reference 

HI 

2009-12-26 

H 

4.2 

3 

162 

33.8 

191' 5 

102.4-255.7 

no 

1 

Paper I 

K1 

2010-11-15 

K s 

10.0 

1 

384 

64.0 

2075 

89.9-263.0 

yes 

0 

this work 

K2 

2012-02-04 

K s 

5.0 

2 

173 

28.8 

571 

105.5-250.4 

yes 

0 

this work 

K3 

2012-11-04 

K s 

0.2 

20 

491 

32.9 

1072 

94.9-260.5 

no 

2 

this work 

K4 

2012-11-06 

K s 

0.2 

20 

231 

15.5 

10/2 

116.3-262.3 

no 

1 

this work 


Notes. Each of the five epochs of observation is assigned an identifier (ID) used throughout this work. The table lists the epoch of observation, 
filter band, detector integration time (DIT) for a single exposure, number of exposures co-added per stored frame, number of usable frames, total 
integration time, full width of the field of view (FoV), range of parallactic angles covered, whether or not the target star was saturated in the science 
images (sat.), the number of PSF reference stars, and the publication describing the data in detail. 


and disks. For point sources, this effect is easy to determine (e.g., 
Lafreniere et al. 2007; Brandt et al. 2013), but the subtraction ef- 
fects are non-trivial for extended sources such as disks and can 
affect the apparent morphology of the source. Forward modeling 
is a powerful tool for interpreting these images and extracting 
the disk geometry (e.g., Thalmann et al. 2011, 2013; Milli et al. 
2012 ). 

One transitional disk that has been the focus of significant 
attention during the past few years is LkCa 15, a Sun-like host 
to a transitional disk with an inner hole the size of our so- 
lar system (~50 AU, e.g., Pietu et al. 2007). Despite the large 
gap, LkCa 15 displays a significant NIR excess (Espaillat et al. 
2007) and residual millimeter emission from small orbital radii 
(Andrews et al. 2011), which implies that an inner AU-sized, 
optically thick disk exists within the gap in addition to the outer 
disk (Espaillat et al. 2008; Mulders et al. 2010). Sparse aperture 
masking observations have revealed an extended structure within 
the gap, which has been interpreted as a possible accreting pro- 
toplanet (Kraus & Ireland 2012). 

Recent simulations explore the structures induced into proto- 
planetary disks by individual embedded planets and predict their 
observable signatures at NIR, mid-infrared, and sub-millimeter 
wavelengths (Zhu et al. 2011; Dodson-Robinson & Salyk 2011; 
Jang-Condell & Turner 2013; de Juan Ovelar et al. 2013). Kraus 
& Ireland (2012) note that their planet candidate cannot be the 
body sculpting the wide gap, suggesting that additional planets 
may reside in the gap. 

In 2010, we reported the first spatially resolved imaging of 
the LkCa 15 gap in the near infrared (Thalmann et al. 2010, here- 
after "Paper I”). The observations confirmed the presence and 
size of the gap as being broadly consistent with that inferred by 
the spectral energy distribution (SED) and seen at longer wave- 
lengths but also indicated a possible offset of the gap center from 
the location of the central star. This implied an eccentric gap 
edge; the likes of which have been suggested as a possible indi- 
cation of planets within the gap, shepherding the disk material 
through their gravitational influence (e.g., Kuchner & Holman 
2003; Quillen 2006). The observations also revealed a strong 
brightness asymmetry between the northern and southern parts 
of the disk, which by itself could either be interpreted as prefer- 
ential forward scattering from optically thin material at the near 
side of the disk or reflection from an optically thick surface at the 
far side of the disk. The disk orientation derived by Pietu et al. 
(2007) on the basis of asymmetries in their millimeter interfer- 
ometry data favors the former scenario, whereas Jang-Condell & 
Turner (2013) assume the latter scenario. 

Overall, the very small angular scale of the disk and the pro- 
hibitive brightness contrast limited the amount of quantitative 
results that could be gleaned from the ADI images in Paper I di- 
rectly, which necessitated the acquisition of deeper high-contrast 


imaging data and a comprehensive forward-modeling effort as a 
next step. 

Here, we present four new epochs of NIR high-contrast 
imaging of the LkCa 15 disk in the K s band. Due to their su- 
perior Strehl ratio, the A^-band data offer cleaner disk images 
than the previously used //-band data. Furthermore, the avail- 
ability of several epochs of observation provides a more robust 
foundation for interpreting the disk morphology through com- 
parison of consistencies and scatter between the epochs. To ex- 
tract quantitative results on the disk geometry from the imaging 
data, we generated an extensive parametric grid of model disks 
as described in Mulders et al. (2010), calculated their appear- 
ance in scattered light, forward-modeled them through the ADI 
process, and compare them to the data using &x 2 metric. 

In the following, we first describe our observations in detail, 
followed by a description of the data reduction procedure. We 
then discuss our modeling efforts to derive the properties of the 
disk and its inner hole, and subsequently show the results of this 
modeling. Finally, we discuss the implications of our results and 
present our conclusions. 

2. Observations 

This work is based on four epochs of //.-band high-contrast 
observations taken with Gemini NIRI (Hodapp et al. 2003) 
from 2010 to 2012 (Gemini Science Programs GN-2010B-Q-93, 
GN-201 1B-Q-36, and GN-2012B-Q-94). In all cases, the adap- 
tive optics was used to deliver diffraction-limited images, and 
the image rotator was operated in pupil-tracking mode to enable 
data reduction with ADI (Marois et al. 2006). The plate scale 
was 07020 per pixel. We hereafter refer to these four epochs 
of A's-band observation as K1-K4. Likewise, we assign the la- 
bel H 1 to the earlier H - band data we published in Paper I, which 
were taken on Subaru HiCIAO (Suzuki et al. 2010) at a plate 
scale of 070095. An overview of all observations and their nu- 
merical parameters is provided in Table 1 . 

The K s -band observations largely share the same observation 
strategy and instrumental setup, which renders them well-suited 
to a systematic analysis. The only change in strategy was the 
adoption of short detector integration times (DIT) in runs K3 
and K4, which allowed for unsaturated imaging of the host star 
LkCa 15 in the science data. For this purpose, the NIRI detec- 
tor was operated in co-add readout mode to avoid large read- 
out overheads. Furthermore, the detector area to be read out was 
windowed from the full 1024 2 pixels down to 256 2 (K2) and 
512 2 pixels (K3, K4), since the astrophysical area of interest for 
this study lies within a radius of ~1" from LkCa 15. 

Due to scheduling constraints and technical downtime, the 
observing runs vary in duration and continuity. The hour angle 
and parallactic angle coverage is illustrated for each run in Fig. 1 . 
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Fig. 1. Parallactic angle coverage of the five observing runs. The num- 
bering of the vertical axis is accurate for epoch HI; all other epochs 
have been translated upwards for visibility. 


While K1 comprises the longest integration time and the fewest 
discontinuities among the A' s -band runs, the target star is satu- 
rated in the science frames, which limits its usefulness for accu- 
rate forward-modeling of the disk (cf. Sect. 4). Furthermore, the 
field rotation stagnates throughout the last hour of Kl; thus, its 
contribution to the ADI data reduction is minimal. Finally, the 
lack of unsaturated stellar PSFs in the science data reduces the 
accuracy of forward-modeling for Kl and K2. Overall, K3 pro- 
vides the best combination of signal-to-noise ratio (S/N) and re- 
liability among our observing runs. We therefore adopt it as the 
benchmark dataset for our further analyses. 

For the observing runs K3 and K4, shorter observations of 
two other target stars were scheduled immediately before and 
after the science observations of LkCa 15 to be used for PSF ref- 
erence. These reference stars, HD 283240 and V1363 Tau, were 
chosen so as to match LkCa 15 closely in color, magnitude, and 
declinations to reproduce the PSF of the science observations as 
well as possible. However, V1363 Tau turned out to be a close bi- 
nary, and thus, it is unsuitable for a PSF reference. Furthermore, 
vibrations in the telescope (Andrew Stephens, p.c.) caused slight 
elongations of the PSF during those epochs, which were less 
prominent for the reference stars than for LkCa 15. As a result, 
we expect the reference PSFs to be of limited use for PSF sub- 
traction of the science data. These elongations should not have a 
significant impact on ADI-based data reductions, though, since 
those techniques use the science observations themselves as 
PSF reference and the position angle of the elongation is sta- 
ble with respect to the pupil. Likewise, our forward-modeling 
analysis described in Sect. 4 takes this effect into account by 
convolving the model disks with the actual PSFs of the science 
data. 


3. Data reduction 

Despite the efforts of the adaptive optics system, the faint scat- 
tered light from the LkCa 15 transitional disk is still over- 
whelmed by the diffraction halo of its host stars in all datasets. 
Differential imaging methods must be employed to remove as 
much of the star’s light as possible and recover the disk flux. 
Some of the disk flux is irretrievably lost in the process as well, 
rendering the reconstruction of the disk's true appearance from 
the resulting images non-trivial (Thalmann et al. 2011, 2013; 
Milli et al. 2012). As discussed in Thalmann et al. (2013), the 
most rigorous way to infer physical disk properties from such 
data is to generate a parametric grid of plausible numerical 
disk models, calculate scattered-light images from those mod- 
els, and subject the theoretical disk models to the exact same 


data reduction process as the science data. The resulting differ- 
ential images can then be compared to those derived from the 
science data to constrain the range of disk parameters consistent 
with observations. 

A number of differential imaging techniques are available 
for use with ADI data of circumstellar disks, offering a trade-off 
between effective suppression of the stellar halo (“aggressive” 
techniques) and conservation of disk flux (“conservative” tech- 
niques). The optimal choice of technique depends on the disk ge- 
ometry, the circumstances of the observation, and the disk prop- 
erties to be measured. 

For the data at hand, we consider four differential imaging 
techniques, which we tested on the K3 dataset. The resulting 
images are presented in Fig. 2. In order of increasing aggres- 
siveness, the methods are as follows: 

PSF reference subtraction based on principle component anal- 
ysis (PCA; Soummer et al. 2012; Amara & Quanz 2012). 
First, we used the Karhunen-Loeve algorithm to decompose 
the dataset of the PSF reference star into a mean PSF plus an 
orthonormal base of principal component images. The im- 
age area within a radius of 6 pixels (O'/ 1 2) was masked out to 
prevent the star’s PSF core to dominate the definition of the 
principal components, as recommended in Amara & Quanz 
(2012). For each frame in the science dataset, we then sub- 
tracted the mean reference PSF and matched the first n prin- 
cipal component templates to the science frame by least ab- 
solute deviation fitting and subtracted them out. The science 
frames were then derotated and co-added to produce the final 
image. We selected n = 2 on the basis of visual inspection of 
the results for n = 1, ... 9. 

Classical ADI as presented in Marois et al. (2006). An estimate 
of the stellar PSF is built by taking the mean of all science 
frames. The PSF is subtracted from all science frames, which 
are then derotated and co-added. We use mean-based frame 
collapse rather than median in order to improve the residual 
noise characteristics, as recommended in Brandt et al. (2013) 
PCA-ADI as employed by the PynPOINT pipeline in Amara 
& Quanz (2012). This technique is essentially identical to 
the PCA-based PSF reference subtraction method described 
above, except that the science dataset itself is used as the ba- 
sis for the principal component templates rather than an ex- 
ternal PSF reference star. Again, we selected n = 2 principal 
components to be subtracted. 

Conservative LOCI as introduced in Paper I and Buenzli 
et al. (2010). The LOCI algorithm (Locally Optimized 
Combination of Images; Lafreniere et al. 2007) is widely 
used in conjunction with ADI for direct imaging of plan- 
ets and brown-dwarf companions (e.g., Marois et al. 2008; 
Thalmann et al. 2009; Lagrange et al. 2010; Carson et al. 
2013). In its most common form, it is too aggressive to im- 
age astrophysical sources beyond the size of an isolated point 
source effectively; however, by choosing a stricter frame se- 
lection criterion (Paper I) or a larger optimization region area 
(Buenzli et al. 2010), the flux loss can be reduced far enough 
to make the method viable for the imaging of circumstellar 
disks. We refer to this use of LOCI as “conservative LOCI.” 
For our K,-band data, we find LOCI to be more aggressive 
than for the //-band data presented in Paper I - possibly 
due to the coarser plate scale - and therefore choose highly 
conservative LOCI parameters: a frame selection criterion of 
Ns = 5 FWHM and an optimization area of Na = 10 000 PSF 
footprints. See Lafreniere et al. 2007 for detailed description 
of these parameters. 
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Fig. 2. Comparison of four high-contrast data reduction methods applied to the K3 dataset of LkCa 15. The top row a)-d) shows the resulting 
images at a fixed linear stretch of ± 1 .6 x 10 -3 times the stellar peak flux. In order of decreasing conservation of disk flux: a) PCA-assisted reference 
PSF subtraction, b) classical ADI, c) PCA-ADI, d) conservative LOCI. The bottom row e)-h) shows the same images after renormalizing each 
concentric annulus around the star by the standard deviation of the pixel values in the annulus at a stretch of ±4cr. The resulting images resemble 
signal-to-noise maps, though the effective noise level is dominated by disk flux and thus overestimated in the inner Nevertheless, these 

images serve to reduce the dynamic range of the high-contrast images and visualize the characteristic crescent of positive disk flux left behind by 
the differential imaging methods. 


All data were Hatfield-corrected and registered using Gaussian 
centroiding of the target star prior to the application of dif- 
ferential imaging. Since all datasets are either unsaturated or 
mildly saturated, the expected registration accuracy is <0.2 pixel 
(4 mas). 

The dynamic range of the resulting images and the varying 
amounts of flux loss among the four reduction methods ren- 
der the images difficult to compare directly. For this purpose, 
we renormalize the images by dividing the pixel values in each 
concentric annulus around the star by the standard deviation of 
those values. The resulting images resemble S/N maps, though 
the “noise” level is dominated by the disk signal, and thus, over- 
estimated out to ~(X'6 (Fig. 2e-h). Nevertheless, they serve to il- 
lustrate the signal content of the differential images qualitatively. 

The images from all four data reduction methods paint 
a consistent picture of a crescent- shaped source of positive 
flux visually consistent in shape, size, and orientation with 
the //-band image published in Paper I, confirming that we 
are indeed detecting scattered light from the surface of the 
LkCa 15 pre-transitional disk. Perhaps surprisingly, the S/N map 
derived from PSF reference subtraction (Fig. 2e) does not 
differ fundamentally from those made with ADI techniques 
(Fig. 2f-h). Although the reference subtraction method avoids 
self-subtraction of the disk, the net positive disk flux still causes 
the template-fitting routine to overestimate and thus oversubtract 
the stellar halo in the science data. Therefore, none of these 
images yield an unbiased view of the LkCa 15 disk, making 
forward-modeling a necessity for quantitative analysis. 

On the basis of Fig. 2, we adopt classical ADI as our differ- 
ential imaging method of choice for the rest of this work. Due to 
the saturation in epochs K1 and K2, PSF reference subtraction 


is ruled out. Of the remaining methods, classical ADI conserves 
the most disk flux, is numerically transparent and linear, and re- 
quires the least computation time. The latter point is relevant be- 
cause the accuracy of a forward-modeling analysis is limited by 
the sampling of the multi-dimensional model parameter space 
and, thus, by the number of models that can be evaluated in a 
reasonable timeframe. 

Figure 3 shows the results of classical ADI applied to all 
four four A/ -band datasets. Overall, the crescent of scattered 
light from the LkCa 15 transitional disk is consistently repro- 
duced among the epochs. In particular, the position and orienta- 
tion of the edge between the disk gap and the bright side of the 
illuminated disk surface can be measured reliably. The ansae of 
the disk gap are difficult to identify due to oversubtraction and 
low S/Ns, and the faint side of the disk surface is not visible. 

4. Modeling 

4.1. Experimental design 

Since the loss of disk flux in ADI image processing is irre- 
versible, the only robust method of extracting astrophysical in- 
formation on the LkCa 15 transitional disk from the data is 
forward modeling. This means generating a large number of 
plausible disk models by means of a radiative transfer code, sim- 
ulating the observable appearance of those disks in scattered- 
light imaging using a ray tracing code, and finally subjecting 
those images to the same ADI image processing that was applied 
to the science data. The family of disk models that yield ADI im- 
ages consistent with those resulting from the science data can 
then be used to derive the best-estimate values and confidence 
intervals for the physical parameters of the disk. 
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Fig. 3. Comparison of the four Gemini NIRI AT s -band observing runs of LkCa 15. All datasets were reduced with classical ADI. As in Fig. 2, 
the upper row shows the resulting images at a linear stretch of ±1.6 x 10 -3 times the stellar peak flux, whereas the bottom row shows radially 
renormalized versions of the same images at a stretch of ±4cr. a) Kl, b) K2, c) K3, d) K4. Since the target star is saturated in runs K1 and K2, 
their flux normalization is approximate. 


Our forward-modeling setup comprises a total of nine inde- 
pendent free parameters: 

- r, the radius of the disk’s transitional gap outer edge, or 
“wall”, represented by the scale factor / := (56 AU)/r; 

- i, the inclination of its orbital plane; 

- g, the Henyey-Greenstein forward-scattering efficiency of its 
dust grains; 

- w, the “roundness” of the gap wall; 

- s, the vertical scale height of the inner disk, with .? = 1 cor- 
responding to the canonical value derived from the SED (see 
Sect. 4.2); 

- o, the orientation of the projected disk’s rotation axis (mea- 
sured from north to east); 

- c, the flux contrast between the disk and the star; 

- ( x , y), the offsets of the disk’s center from the star parallel 
and perpendicular to the line of nodes, respectively. The line 
of nodes is the intersection of the inclined disk plane with 
the “sky plane” perpendicular to the line of sight; it defines 
the unforeshortened “major axis” of the projected disk. The 
parameter y includes foreshortening; its observed value is 
cos(i) times the physical offset along the disk plane. 

The radiative transfer code does not support eccentric disks. 
However, since we discovered structure indicative of eccentric- 
ity in the LkCa 15 disk in Paper I, it is scientifically interesting 
to test this hypothesis in the analysis at hand. We therefore ap- 
proximate eccentric disks by calculating scattered-light images 
of azimuthally symmetric disks, translating them by a small off- 
set ( x , y) with respect to the star, and rescaling the brightness to 
take into account the new center of illumination. This is done in 
the third stage of our analysis (ADI processing). 

Figure 4 provides a graphical visualization of the LkCa 15 
system architecture as described by our model. Note that the 
sketches are not to scale - the SED predicts the inner truncation 


radius of the outer disk to be of order 50 AU, whereas the inner 
disk is contained at sub-AU radii. As a result, the inner disk is 
inaccessible to our direct imaging efforts and can only impact 
our observations through the shadow it may cast on the wall of 
the outer disk. 

The radiative transfer code used to generate self-consistent 
physical disk models is described in Sect. 4.2, the raytracing 
code employed to generate scattered-light images from the disk 
models in Sect. 4.3, and the ADI processing and^ 2 evaluation 
in Sect. 4.4. 

4.2. Radiative transfer code 

The radiative transfer code used in this paper is MCMax (Min 
et al. 2009), a disk modeling tool that performs 3D dust ra- 
diative transfer in a 2D axisymmetric geometry. The code in- 
cludes full anisotropic scattering and polarization (Min et al. 
2012; Mulders et al. 2013a), making it an ideal tool for inter- 
preting high-contrast images of protoplanetary disks. 

Besides radiative transfer, it solves for the vertical structure 
of the disk using hydrostatic equilibrium and dust settling, which 
yields a self-consistent density and temperature structure. The 
dust scale height is in this case controlled by a scale factor Y, 
which represents the reduction in scale height of the dust com- 
pared to that of the gas. 

The model employed here is based on the optically thick 
inner disk model described in Mulders et al. (2010, hereafter 
M10). The main modifications include: 

Anisotropic scattering. The M10 model did not include scatter- 
ing in the radiative transfer step of the simulation. 

Dust properties. We changed the dust size and composition to 
roughly reflect the observed brightness asymmetry, such that 
the imaging constraints do not significantly alter the SED 
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(a) Front view inner disk 
(not to scale) 


outer disk 
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nodes 
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Jl 
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w= 0 
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Fig. 4. Sketch of the LkCa 15 system architec- 
ture as described by our model. Most of the 
bulk of the outer disk (shown in light grey) re- 
mains unseen; only the starlight reflected from 
its surface (shown in red) is observed with di- 
rect imaging. Anisotropic scattering behavior 
greatly enhances the forward-scattered flux on 
the disk’s near side, creating the bright cres- 
cent in the ADI images. The treatment of wall 
shape in our model is visualized in more de- 
tail in Fig. 5. By default, the outer disk scale 
height h is a fixed rather than a free parameter; 
see Sect. 5.9 and Appendix C. 


fit. The dust properties are described in the next section and 
shown in Table 2. 

Round wall. Steep gradients in surface density are not a natu- 
ral outcome of planet-disk interactions (Lubow & D’Angelo 
2006; Crida & Morbidelli 2007). Rather than a sudden in- 
crease at r, the surface density I gradually increases with 
radius up to a radius R exp after which it follows the surface 
density of the outer disk (E 0 uter R~ l ). We use the exponential 
function described in Mulders et al. (2013b), which is char- 
acterized by a dimensionless width w to described the spatial 
extent of this transition: 


Y.^R < R e x p) — Pouter R CXp I 


1 - R/R { 


exp 


(i) 


This smoothly increasing surface density gives rises to a 
rounded-off wall in spatially resolved images, as shown in 


Fig. 5. The wall extends further in than the anchor point for 
the exponential turnover (R exp ~ [1...3]r), depending on wall 
roundness). Instead, we use the radial peak in the intensity to 
trace the wall location r. Thus, w represents the characteris- 
tic length scale of the exponential decay of the disk’s surface 
brightness at the gap edge. It is a dimensionless parameter 
normalized to the gap radius r. Both the wall location r and 
its roundness w are free parameters. 

Stellar properties. We use the updated stellar parameters from 
Andrews et al. (2011), as listed in Table 2. 

The dust contains two components with different size distribu- 
tions (flsmaii and abig) and different degrees of settling (.? for 
the inner disk and T' sma ii and T'big for the outer disk). The pa- 
rameters T'smaii and 'Fbig are constrained by fitting the mid- 
infrared SED. The variable .v determines the size of the shadow 
cast on the outer disk wall. Although fitting the NIR excess 
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Table 2. Disk parameters for the geometrical model. 


Parameter 

Value 

Teff [K] 

4370 

L, [Lo] 

1.22 

M, [M g ] 

1.01 

d [pci 

140 

R m [AU1 

0.07 

Rou, [AU] 

300+ 

Rgap [AU] 

4.. .56 

A/dust \M q ] 

0.0 1 

A/iust, inner [A/©] 

5 x 10" 11 

p 

V 

composition 

SolaU 

^smaii [/^rn] 

0.1. ..1.5 

^big [/^FTll 

100.. .500 

/big, outer L%J 

99 

YU 

1 small 

0.65 

** 

0.2 


Notes. Fixed parameters are indicated by a dagger (t). The surface den- 
sity profile is defined as £(r) cc r~ p and scaled to the total disk mass. The 
location of the gap edge is approximate due to its smooth nature; dis- 
played here is the peak of the radial intensity (Fig. 5). Opacities are cal- 
culated assuming a grain size distribution f(a) oc a -3 - 5 . Solar composi- 
tion: 12% MgFeSi0 4 , 12% MgFeSi 2 0 6 , 12% Mg 2 Si0 4 , 12% MgSi 2 0 6 , 
15% FeS, and 40% C. Optical constants are from: silicates (Dorschner 
et al. 1995; Henning & Stognienko 1996; Mutschke et al. 1998), carbon 
(Preibisch et al. 1993), troilite (Begemann et al. 1994). 



R [AU] 


Fig. 5. Azimuthally averaged surface brightness profile (top) and corre- 
sponding surface density profile (bottom), demonstrating the effect of 
rounding off the disk wall (w > 0). The solid line shows the best-fit 
model with a wall shape of w = 0.30, as compared to a more vertical 
wall with w = 0.05 (dashed line). The dotted line indicates the “wall 
location” r , corresponding to the radial peak in intensity. 

yields s - 1 (M10), we treat this parameter as free because the 
inner disk scale height is known to vary (Espaillat et al. 2011) 
and may be different for each observing epoch. 

Since we updated both the stellar parameters and the dust 
composition since the M10, it was needed to refit all parameters 
to the observed SED (Fig. 6.). All fit parameters are shown in 
Table 2. Note that parameters s, w, and g are also used in the 
radiative transfer step but are free parameters. 

4.2.1. Grain size 

The absorption and scattering properties of the dust are set by the 
grain size, structure, and composition (e.g.. Van der Hulst 1957). 



Fig. 6 . Spectral energy distribution of LkCal5. Displayed with the 
solid line is the best-fit disk model described in Table 2. References: 
(UBVRI) Kenyon & Hartmann (1995); (JHK) Skrutskie et al. (2006); 
Spitzer / IRAC (Luhman et al. 2010); (IRAS) Weaver & Jones (1992); 
(sub-mm) Andrews & Williams (2005); (mm) Pietu et al. (2006). 

The main diagnostics of these properties in unpolarized images 
are the wavelength-dependent albedo (Fukagawa et al. 2010; 
Mulders et al. 2013a) and phase function (Duchene et al. 2004). 
Because the latter is a highly non-linear function of grain size, 
we use the Henyey-Greenstein phase function with assymetry 
parameter g (Henyey & Greenstein 1941) to prescribe the phase 
function of our particles rather than varying the grain size. In do- 
ing so, we avoid having to refit the SED for each different grain 
size, as the grain properties can be kept fixed while varying g. 

The grain composition is based on a condensation sequence 
for solar system elemental abundances as described in Min et al. 
(2011); see Table 2. We assume the grains have an irregular 
shape, parametrized by a distribution of hollow spheres (DHS) 
with a vacuum fraction ranging from 0 to 0.7 (Min et al. 2008). 
With an MRN size distribution (Mathis et al. 1977) from 0.1 
to 1.5 micron, the albedo of these grains is ~0.5 in IGband with 
a phase function close to g = 0.5. 

4.3. Scattered-light image simulation 

The scattered light images are calculated by integrating the for- 
mal solution to the radiative transfer equations, using the den- 
sity and temperature structure from the Monte Carlo simulation 
and the dust opacities as described above. The local scattered 
field is calculated in 3D using a Monte Carlo approach (see Min 
et al. 2012 for details). This scattered field includes a contribu- 
tion both from the star 1 and from the disk, which can be signifi- 
cant at NIR wavelengths. 

First, the images are calculated for a given inclination angle i 
in spherical coordinates, preserving the radial grid refinement 
at the gap outer edge. Finally, this image is mapped onto a 
Cartesian grid matching the pixel scale of the observations, 
yielding a scattered light image that can be further processed 
by ADI. By employing a scale factor / for the field of view, 
we generate multiple Cartesian images with the same resolution 
but a different field of view from the same spherical coordinates 
image. Hence, a larger field of view (large /) corresponds to a 
smaller gap. This greatly decreasing runtime of the entire grid 
by skipping the radiative transfer and raytracing step that would 
otherwise be associated with changing the gap radius R. 

1 The star is treated as a disk of uniform brightness with radius 

R= 1.65 R e . 
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4.4. Forward-modeling of ADI observations 

4.4.1 . Theory 

The raw science data I used as input for the ADI process can be 
described as a sum of the stellar PSFs S and the scattered-light 
images of the disk D\ 

i = s + n. ( 2 ) 

Since classical ADI using mean-based frame combination is lin- 
ear (Marois et al. 2006), the resulting output image O produced 
by ADI processing of the imaging data I is equal to the sum 
of the ADI-processed stellar PSFs and the ADI-processed disk 
images: 

O = ADI(J) = A DIOS’ + £>) = A DIOS) + ADI(£>). (3) 

Since the stellar PSF remains largely static throughout the sci- 
ence data, ADI effectively removes the bulk of the stellar flux, 
leaving behind a halo of residual speckle noise. Given good 
adaptive optics performance (stable Strehl ratio) and enough 
field rotation (ideally across several resolution elements at all 
considered radii), the noise is well behaved and approximately 
Gaussian in concentric annuli around the star. That is, 

O = ADI(D) + noise. (4) 

Thus, applying the ADI process to a noise-free model disk im- 
age lD m od yields an output image ADI(D mo( j) that can be directly 
compared to the resulting image from the science data. In the 
ideal case, subtracting the ADI-processed model image from the 
science data should leave behind only pure noise: 

<3-ADI(D mo d) = ADI(D) + noise - ADI (D mo d) 

= [ADI(D) - ADI (D m od)] + noise 

I 

= noise. (5) 

Finally, assuming that the noise behaves in a Gaussian manner, 
one can define a x 2 metric to measure the goodness of fit for a 
given disk model: 

a v [0-ADiaD m0 d)] 2 

*mcd = 2, ^2 < 6 > 

where cr is a map of the local standard deviation of the noise. 

4.4.2. Noise estimation 

For highly inclined and well resolved circumstellar disks, the ra- 
dial noise profile can be measured in sectors of the final science 
image unaffected by the disk flux (e.g., Thalmann et al. 2011, 
2013). However, the crescent of scattered light from the LkCa 15 
disk, as well as the surrounding oversubtraction regions, dom- 
inate our ADI images at all position angles, which renders it 
impossible to measure an unbiased noise profile. During the 
exploratory phase of our experiment, we therefore started out 
with an estimated noise profile. As soon as we had located 
well fitting models for each epoch, we iteratively redefined the 
noise profile for each dataset as the standard deviation of the 
pixel values in concentric annuli measured in the residual image 
O - ADI(Dbest) for the best-fitting model disk ©best- Such an 
a posteriori definition of the noise profile yields a reduced x 2 
of ~1 for the best-fit model by design and, therefore, cannot 
be used as an unbiased measure of absolute goodness of fit. 
However, it provides a useful measure of relative goodness of 


fit, allowing us to converge on a best-fit set of model parame- 
ters and to define confidence intervals around those parameters 
through ax 2 threshold. 

The noise profiles obtained with this method are furthermore 
used to generate more meaningful S/N maps than the ones shown 
in Figs. 2 and 3, since the disk flux no longer dominates the def- 
inition of the noise profile. The S/N map of the residual image, 
[O - ADI(23best)]/0% provides a useful visual representation of 
the absolute goodness of fit, as well as a “sanity check” on the 
assumed a posteriori noise map. Ideally, it should look like ran- 
dom noise with no coherent structure from the disk left behind. 
Similarly, the S/N map of the unsubtracted science ADI image, 
O/cr, provides a measure of the significance of the retrieved 
disk signal and, thus, of the scientific information content of the 
dataset. 

4.4.3. Implementation 

We perform the forward-modeling of ADI observations for our 
model disk images with a custom IDL code, which comprises 
the following steps: 

- Load a simulated scattered-light disk model for a given set 
of parameters (r, i, g, w, s). 

- If ( x , y ) are not zero, adjust the brightness distribution of the 
image to represent illumination by the off-centered star. For 
this purpose, we multiply each pixel by the square of its dis- 
tance from the image center, then divide by the square of the 
distance from the new, offset position of the star. When cal- 
culating distances, we assume that all disk flux originates in 
the midplane; i.e., the projection effect from the plane’s in- 
clination is taken into account but not the finite scale height 
of the scattering disk surface. 

- For each exposure n in the science data, generate a corre- 
sponding “exposure” of the model disk image. The position 
angle of the model disk in each frame is the input parame- 
ter o plus the parallactic angle of the corresponding science 
exposure. As part of the image rotation, the intended posi- 
tion of the star ( x , y ) in the model image is mapped onto the 
center of the model exposure. This way, the model images 
need only be resampled once, minimizing aliasing effects. 

- Convolve each model exposure with a 15 x 15 pixel image 
(~5 x 5 resolution elements) of the unsaturated PSF core of 
the corresponding science exposure. This gives the simulated 
disk data the appropriate spatial resolution while avoiding 
contamination of the model disk with the physical disk struc- 
ture present in the science exposure at larger separations. In 
the case of saturated science data (epochs Kl, K2), we sub- 
stitute a single external PSF for all exposures. 

- Perform classical ADI on the set of model exposures, yield- 
ing a noise-free ADI output image for the model disk. 

- Bin down both the model output image and the science 
output image by a factor of 3 in both dimensions for the 
purpose of x 2 evaluation. Since the diffraction-limited res- 
olution A/D is 57 mas = 2.8 pixels, neighboring pixels are 
always strongly correlated. After the binning, each pixel in 
the binned image can be treated as an independent data point. 
In practice, however, pixels may still be correlated to a cer- 
tain degree due to large-scale structures remaining in the 
residual image. We retain the unbinned images for visual- 
ization purposes. 

- Match the model output image to the science output im- 
age with a least-squares optimization of the model output 
image’s intensity scale factor. We make use of the linfit 
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function in IDL, which is weighted with the a posteriori 
noise map cr. This defines the final free parameter of the disk 
model, the disk/star brightness contrast c. 

- We subtract the model output image from the science out- 
put image and calculate the^ 2 value of the residual image as 
per Eq. (6), using the a posteriori noise map cr. We restrict 
the evaluation region to the annulus between the radii O'.' 2 
and O'.' 8, thus excluding the center dominated by the stel- 
lar PSF core. The evaluation region comprises N = 516 
binned pixels, each of which roughly corresponds to a res- 
olution element. 

4.4.4. Confidence intervals 

The best-fit disk model is defined as the set of input parameters 
yielding the minimum y 2 value for each epoch of observation. 
We determine confidence intervals around those parameters by 
characterizing the family of “well fitting” solutions bounded by 
a threshold of^ 2 < y^ jn + Ay 2 . As in the case of our analysis 
of the HIP 79977 debris disk (Thalmann et al. 2013), we find 
that the pixels in our residual images do not have fully indepen- 
dent normal errors, as evidenced by traces of large-scale struc- 
ture that remain visible in the S/N maps. Applying the canon- 
ical y 2 thresholds would not take into account those remaining 
systematic errors and, thus, overestimate the confidence in the 
best-fit parameters. We therefore choose a more conservative 
threshold of A^ 2 = 32 ~ V2 N, which corresponds to a lcr de- 
viation of the^ 2 distribution with N = 516 degrees of freedom. 
In other words, this threshold delimits the family of solutions 
whose residual images are consistent with the a posteriori noise 
map derived from the best-fit residual at the 1 cr level. 

Assuming that the disk image does not change significantly 
with time, a global best-fit solution and well fitting solution fam- 
ily can be obtained by co-adding the four^ 2 maps and adjusting 
the threshold A^ 2 by a factor of V4 = 2. The assumption is jus- 
tified for most model parameters, since the disk gap is known 
to lie at a separation of ~50 AU (Espaillat et al. 2008; Andrews 
et al. 2011), which implies an orbital time scale of several hun- 
dred years. Thus, we do not expect large-scale disk properties 
like the inclination, gap radius, or eccentricity of the outer disk to 
evolve measurably between our observing epochs. On the other 
hand, the inner component of the LkCa 15 pre-transitional disk 
orbits at sub-AU separations and may therefore evolve at a time 
scale of months. The model parameter s, which describes the 
vertical scale height of the inner disk and thus the width of the 
shadow band on the outer disk wall, must therefore be consid- 
ered variable between Kl, K2, and the pair of almost contem- 
porary epochs (K3, K4). In any case, we present our results for 
both individual epochs and the combined analysis to allow di- 
rect comparison. As reported in Sect. 5, no significant temporal 
variation is observed in any of the model parameters, including 
the inner disk scale height This validates the use of the global 
best-fit solution derived from all four epochs. 

5. Results 

5.1. Overview 

In the early stages of this analysis, we explored parameter grids 
involving two or three model parameters at a time, keeping the 
other parameters fixed, and iterated this process until a stable 
y 2 minimum was located for each observing epoch. These pre- 
liminary best-fit (PBF) solutions are listed in Appendix A. While 
these models yielded visually convincing residuals and provided 


a necessary starting point for further analyses, they lacked con- 
fidence intervals and were not proven to be the global minima 
of the y 2 landscape. We do not document the specifics of these 
early searches in this work, since they are rendered obsolete by 
the further analyses presented in the following sections. 

The most robust way to establish confidence intervals 
is a comprehensive, finely-sampled exploration of the nine- 
dimensional model parameter space. However, the radiative 
transfer simulation and the raytracing are computation-intensive 
procedures, which renders such a brute-force approach imprac- 
ticable. We therefore seek to reduce the dimensionality of the 
parameter space as much as possible beforehand. 

Simply exploring each parameter individually in one- 
dimensional analyses around the best-fit solution is not a valid 
approach, since some parameters are strongly covariant with 
each other. For instance, the most striking feature of the ADI im- 
ages of LkCa 15 is the sharp transition from the dark disk gap 
to the brightest edge of the reflected light on the near side. For 
any given value of the gap radius r. there is only a narrow range 
of inclinations i that project the gap edge onto the correct appar- 
ent separation from the star to match the data. However, a good 
match can be achieved over a much wider range of inclinations 
if the radius is adjusted proportionally. The well fitting solution 
family therefore lies along a diagonal in the (r, i) plane and is 
ill-described by its cross-sections with the r and i axes. 

Nevertheless, we find that three parameters can be safely de- 
coupled from the full parameter space and that one additional 
parameter requires only minimal sampling. This leaves only a 
manageable five-dimensional parameter space to be explored 
thoroughly by brute-force calculations. 

- c: the flux scaling of the model disk is set by least-squares 
optimization using the linfit function of IDL as part of 
the evaluation code (cf. Sect. 4.4.3). Thus, the disk/star con- 
trast c does not contribute a dimension to the parameter grid 
to be sampled. 

- x , o: the gap center offset x along the line of nodes (the “ma- 
jor” axis of the projected gap ellipse) and the orientation an- 
gle o are the only two parameters whose effects on the disk 
image are not laterally symmetric. Thus, we expect them to 
have well defined optima largely independent of the other 
model parameters. 

- s: our preliminary analyses surprisingly indicated that the 
achievable y 2 monotonously decreases with decreasing in- 
ner disk dust scale height reaching its minimum at s = 0. 
The value s = 1 corresponds to the hydrostatic equilibrium 
value of the inner disk scale height, which is supported by 
observations (Espaillat et al. 2008, 2010). Fractional values 
of s <K 1 correspond to a razor-thin disk, which is not physi- 
cally plausible. The special case of s = 0, on the other hand, 
matches another physically plausible scenario: either the in- 
ner disk is inclined with respect to the outer disk, such that 
most of the inner disk’s shadow bypasses the outer disk, or 
the inner disk is radially optically thin as described in M10; 
in both cases, leaving the outer disk fully illuminated. We 
therefore limit the sampling in our analysis to the discrete 
values s = {0, 0.5, 1} to verify this result. 

5.2. Constraints on the orientation angle o and on the center 
offset along the line of nodes x 

To determine the best-fit values and confidence intervals for x 
and o, we run three-dimensional parameter grids by vary- 
ing (x, o, r ) and keeping all other parameters fixed at their 
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Fig. 7. Constraints on disk center offset x along the line of nodes. The 
plot shows the excess of the best-fit^ 2 for a given x with respect to the 
minimal in achieved, which is normalized by the threshold value A^ 2 . 
The color-coded curves represent each individual epoch, whereas the 
thick black curve is obtained by combining all four epochs. The dashed 
line marks the case of a laterally symmetric disk ( x = 0), whereas the 
dotted line marks the tentative offset postulated in Paper I. Positive val- 
ues of x represent offsets in the direction of the western ansa. The bot- 
tom panel shows the well fitting range and the best-fit value for all five 
analyses. 


preliminary best-fit values. The disk radius r is included a grid 
parameter to allow the positions of the two ansae along the line 
of nodes (±r + x) to be optimized independently of each other. 

Figures 7 and 8 show the resulting best-fit values and do- 
mains of the well fitting solution families for x and o, respec- 
tively. The families from the four epochs appear as well-behaved 
X 2 curves with unique minima. While the o values from different 
epochs agree well, the x values exhibit some scatter, in particular 
in the lower-quality epochs K2 and K4. 

We find that the best-fit offsets x consistently lie on the pos- 
itive (western) side of the star in all four epochs, though the 
offsets for the two high-quality datasets K1 and K3 are consis- 
tent with zero. Adding up the four^ 2 maps yields global best- 
fit values and well fitting intervals of x = +24 [+7, +45] mas 
(corresponding to +3 [+ 1 , +6] AU in physical distance) and 
o = 150 [147. 153] degrees. A summary of all numerical results 
from this analysis is included in Table 3. 

In Paper I, we postulated a tentative disk center offset of 
x = +64 mas with a phenomenological error bar of ±6 mas 
on the basis of the HI dataset. This error did not take into ac- 
count the bias and systematic uncertainties from the ADI re- 
duction and therefore overestimated the accuracy of the mea- 
surement. However, the measured value is consistent with those 
of our lower-quality datasets K2 and K4, whose S/N maps ex- 
hibit a similarly "patchy” structure like that of HI. The spread 
in x values may therefore be due to a bias caused by marginal 
lateral constraints on the gap; i.e., the optimization of x might 
be dominated by large-scale residual noise rather than by the 
actual disk flux in the lower-quality datasets. The fact that all 
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Fig. 8. Constraints on the orientation angle o of the disk’s rotation axis. 
The plot shows the excess of the best-fit x 1 f° r a given o with respect 
to the minimal y 2 . achieved, which is normalized by the threshold 
value tsx 1 ■ The color-coded curves represent each individual epoch, 
whereas the thick black curve is obtained by combining all four epochs. 
The angle o is measured counter-clockwise from north to east. The bot- 
tom panel shows the well fitting range and the best-fit value for all five 
analyses. 


five measurements of x have the same sign suggests an underly- 
ing physical asymmetry that may be exaggerated by ADI under 
low S/N conditions. 

Thus, our analysis overall tentatively confirms the existence 
of a positive disk center offset in x direction, though at a signifi- 
cantly smaller magnitude than proposed in Paper I. 

Note that the bright crescent apparent in the ADI images is 
generated by the near side of the disk for all our best-fit mod- 
els, where the scattered light is greatly enhanced by anisotropic 
forward scattering (“near-side scenario”). This agrees with the 
predictions of Pietu et al. (2007) based on asymmetries in their 
millimeter interferometry images. Furthermore, our orientation 
angle agrees very well with their value of o = 150.7° ± 0.4°. 

In Paper I, we demonstrated that it is possible to make the far 
side of the disk rim brighter than the near side within the param- 
eter space of our model, allowing for the “far-side scenario” as 
an alternative explanation for the bright crescent in the ADI data. 
In this scenario, the scattering anisotropy is minimal ( g ~ 0), and 
the contrast is instead generated by self-shadowing of the near- 
side gap wall, which leaves only the illuminated far-side wall 
exposed to imaging. To explore this possibility, we have run a 
limited parameter grid of models with o = 230°, g = 0, coarse 
sampling of s and w, and fine sampling of i, /, and y. The latter 
three parameters are most likely to influence the near-to-far side 
contrast. 

The results are shown in Fig. 9. While the far-side scenario 
is capable of achieving a contrast comparable to the one in the 
near-side scenario, it requires the visible scattered light to be 
tightly constrained to the gap wall. The resulting image is a 
poor match for the observed radially extended disk morphol- 
ogy. The best^ 2 attained with the far-side model exceeds that of 
the near-side model by 18A^ 2 . We therefore conclude that the 
bright crescent of scattered light in the LkCa 15 disk represents 
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Fig. 9. Comparison of the best-fit solutions for the near-side (g = 0.7, 
o = 150°) and the far-side (g = 0, o = 230°) scenarios for the 
K3 dataset, a) Image of the unconvolved best-fit model near-side disk 
model at logarithmic stretch, b) Same for the best-fit far-side model. 
Note that the dark lane on the faint side is not generated by the shadow 
cast from the inner disk in this configuration, but rather by the opti- 
cally thick bulk of the outer disk blocking the view onto the illuminated 
face of the gap wall. The two parallel arcs of light are scattered light 
spilling around the edges of the near-side wall, c) S/N map of the resid- 
ual image after subtracting the ADI-processed near-side model image 
from the ADI-processed data at a linear stretch of ±5cr. d) Same for 
the far-side model, e) Difference of the images d)-c), at a linear stretch 
of ±2cr. The^ of the far-side model exceeds that of the near-side model 
by 18A_^ 2 . We therefore conclude that the bright crescent in the data 
represents the strongly forward-scattering near side of the disk. 


starlight forward-scattered off the near-side surface of the outer 
disk. 


5.3. Brute-force optimization and best-fit solutions 

Having obtained best-fit solutions for x and o, we perform a 
brute-force exploration of the remaining parameter space with 
fine sampling in /, i, g , w , y, discrete sampling in s, and least- 
squares fitting in c. Table 3 defines the range and step size of the 
sampling in all parameter dimensions. 

We find well defined^ 2 minima for all four epochs. Table 3 
lists the model parameters corresponding to these best-fit solu- 
tions and the parameter ranges of the well fitting solution family 
obtained with a threshold of x 2 + A^ 2 as defined in Sect. 4.4.4. 
A better characterization of the well fitting solution family and 
the covariances between the model parameters is provided by 
the^ 2 contour plots in the following subsections. 

Figure 10 presents a visual overview of the best-fit model 
solutions and their agreement with the data for each epoch of 
observation. The first row displays the classical ADI images de- 
rived from the data, while the second row shows the classical 
ADI images generated from the best-fit model disk images with 
the same set of parallactic angles. The third row represents the 
subtraction residuals of the first two rows, demonstrating that the 
overall morphology of the data is well explained by the models. 
The residuals are largely unstructured and well behaved; a coher- 
ent arc reminiscent of the disk structure is retained only in K2. 
There are no consistently recurrent features apparent in the four 
residual images that would suggest observable disk structure be- 
yond the scope of the model parameter space, such as spirals, 
clumps, gaps, or planets. 

The fourth row exhibits the unprocessed simulated scattered- 
light images of the best-fit model disks in logarithmic stretch. 
The disk architectures are roughly consistent between the 
epochs. Note that the position of the disk center is consistently 
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Fig. 10. Best-fit models for the LkCa 15 disk. The columns represent the four epochs of observation, K1-K4. a)-d) ADI output images from 
science data; e)-h) ADI output images from the simulated model disks; i)-l) science-model subtraction residuals at the same linear stretch; 
m)-p) unprocessed model disk images at logarithmic stretch; q)-t) S/N maps of the science ADI output images at a linear stretch of [-5cr, 5crl; 
u)-x) S/N maps of the subtraction residuals at the same stretch. 
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Table 4. Information content of the four datasets, as measured by the 
reduced^ 2 of the a posteriori S/N map of the raw ADI disk image. 


K1 

K2 

K3 

K4 

Reduced^ 2 3.9 

2.1 

5.4 

2.5 


Notes. A higher reduced^ 2 implies an ADI image that is more inconsis- 
tent with being pure noise, or that has more disk flux measurable above 
the noise threshold. 

offset from the star, bringing the bright near side of the disk rim 
closer to the star. In our model parametrization, this corresponds 
to a negative value of y. This observation is further discussed 
below in Sect. 5.9. 

The fifth and sixth rows of Fig. 10 show the S/N maps cor- 
responding to the original ADI output images (first row) and 
subtraction residuals (third row), using the a posteriori radial 
noise map derived from the residuals (third row), as discussed in 
Sect. 4.4.2. This reduces the dynamic range of the images, allow- 
ing residual structures at various radii to be compared directly. 
The images of the fifth row further offer an estimate of the over- 
all information content of each observing run. Since the original 
disk flux can be assumed as approximately constant, the S/N at 
which the disk is detected at a given location is a measure of 
sensitivity and effectiveness of noise suppression. 

We can obtain a numerical measure of the information con- 
tent in each dataset by calculating the reduced x 2 of the raw 
ADI images, which is defined as the RMS of the S/N map di- 
vided by the number of pixels in the evaluation area. The result- 
ing numbers are given in Table 4. This confirms our choice of K3 
as the most reliable dataset and explains the inferior performance 
of K2 in Fig. 10. 

5.4. Constraints on the wait shape parameter w 

Figure 1 1 presents the minimum x 2 achieved in our brute-force 
grid search as a function of the gap wall shape parameter w. The 
best-fit solutions of all four epochs lie within the narrow interval 
of w = [0.27, 0.31]. While two individual epochs cannot exclude 
the possibility of an abrupt “vertical” wall ( w ~ 0), the combined 
analysis yields a well fitting range of w = [0.19,0.36], which 
clearly favors an extended, “fuzzy" radial transition from gap to 
outer disk. 

Figure 1 2 illustrates the difference between a round wall (s = 
0.30) and a vertical wall (5 = 0.05) for the epoch K3. We find 
no significant codependence of the best-fit value of w with other 
model parameters. 

5.5. Constraints on inner disk settling s 

Figure 13 shows a plot of the minimal^ 2 achieved in our brute- 
force grid search for each of the three considered values of the 
inner disk settling parameter s = jO, 0.5, 1] compared to the 
global x 2 minimum. As expected on the basis of our prelimi- 
nary analyses (cf. Sect. 5.1), the global best fit is achieved with 
s = 0 for all epochs, whereas the best fits attainable with s = 1 
exceed the global^ 2 minimum by 3-7 Ax 2 (TO A^ 2 for all four 
epochs combined). 

The figure suggests that the criterion for a well fitting so- 
lution (x 2 < X^ n + At 2 ) i s only achievable with s <K 1. 
Taken at face value, this implies an unphysically small scale 
height for the inner disk, given that 5=1 represents the ex- 
pected hydrostatic equilibrium value that is compatible with the 



Wall shape parameter w 

Fig. 11. Constraints on the wall shape parameter w. The plot shows the 
excess of the best-fit x 2 for a given w with respect to the minimal x^ n 
achieved, which is normalized by the threshold value A^ 2 . The color- 
coded curves represent each individual epoch, whereas the thick black 
curve is obtained by combining all four epochs. The bottom panel shows 
the well fitting range and the best-fit value for all five analyses. The 
global best fit is achieved for w = 0.31, suggesting that the gap wall 
of the outer disk is “fuzzy”, extending over a significant range of radii 
rather than being sharply constrained to a given radius (w ~ 0). 



Fig. 12. Comparison of the best-fit solutions for w = 0.30 and w = 0.05 
for the K3 dataset, a) Image of the unconvolved best-fit model disk with 
tn = 0.30 at logarithmic stretch, b) Same for w = 0.05. c) S/N map of 
the residual image after subtracting the ADI-processed best-fit model 
image with w = 0.30 from the ADI-processed data at a linear stretch 
of +5cr. d) Same for w = 0.05. e) Difference of the images d)-c) at a 
linear stretch of ±2o\ The^- 2 of the w = 0.05 model exceeds that of the 
tn = 0.30 model by 0.85 A^ 2 . 


observed near-infared excess. However, an inner disk inclined 
with respect to the outer disk would leave most of the outer 
disk wall unshadowed and, thus, provide a physical scenario for 
5 = 0. We note that the alternative scenario of an optically thin 
dust halo in place of the inner disk, as proposed in Espaillat 
et al. (2007) and investigated in Mulders et al. (2010), is also 
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Inner disk settling parameter s 

Fig. 13. Constraints on the inner disk vertical settling parameter s. The 
plot shows the excess of the best-fit for a given i with respect to the 
minimal ;^ 2 in achieved, which is normalized by the threshold value . 
The color-coded curves represent each individual epoch, whereas the 
thick black curve is obtained by combining all four epochs. In all cases, 
the best fit is achieved with s = 0, or with a completely absent inner 
disk shadow. 



Fig. 14. Comparison of the best-fit solutions for j = 0 and s = 1 for the 
K3 dataset, a) Image of the unconvolved best-fit model disk with s = 0 
at logarithmic stretch, b) Same for s = 1. Note the lane of shadowing 
from the inner disk visible on the far-side gap wall, c) S/N map of the 
residual image after subtracting the ADI-processed best-fit model image 
with s = 0 from the ADI-processed data at a linear stretch of ±5cr. 
d) Same for s = 1. e) Difference of the images d)-c) at a linear stretch 
of ±2o\ The^ 2 of the s = 1 model exceeds that of the r = 0 model by 
7.4 A* 2 . 
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Fig. 15. Constraints on the outer disk wall radius r, as represented in our 
model by the scale factor / = 56 AU/r. The plot shows the excess of the 
best-fit^ 2 for a given r with respect to the global minimum which 
is normalized by the threshold value A^ 2 . The color-coded curves rep- 
resent each individual epoch, whereas the thick black curve is obtained 
by combining all four epochs. The bottom panel shows the well fitting 
range and the best-fit value for all five analyses. 


cannot be distinguished in our data. However, note that the effect 
of shadowing from the inner disk in the s = 1 case is not re- 
stricted to the faint far-side gap wall, as one might expect. While 
the near-side wall would be self-shadowed in the case of a verti- 
cal cutoff (w = 0), its upper half remains exposed to direct view 
in the case of a tapered disk edge ( w > 0), which is favored by 
our analysis (cf. Sect. 5.4). Thus, the inner disk shadow visibly 
truncates the inner edge of the tapered disk wall, hardening the 
gradient between the illuminated disk surface and the gap at all 
azimuths. As a result, our analysis is more sensitive to the degree 
of shadowing than the low S/N of the far-side gap edge would 
suggest. 

Unlike all other model parameters, s pertains to the inner 
rather than the outer disk, which has an orbital time scale of 
months rather than decades. As a result, it is the only parame- 
ter for which one may expect astrophysical variation from one 
observing epoch to the next. However, we find consistent results 
for s = 0 in all datasets. 

Within the parameter space of our modeling effort, we there- 
fore conclude that an inner disk inclined with respect to the outer 
disk is the most likely scenario for LkCa 15. 


compatible with s = 0, but the existence of such a halo is diffi- 
cult to justify from a physical point of view (Espaillat et al. 2010; 
however, see Krijt & Dominik 2011). 

Figure 14 illustrates the differences between the best-fit so- 
lutions for s = 0 and s = 1. The shadow of the inner disk on 
the outer disk wall is visible as a dark ribbon in the 5 = 1 model 
image (Fig. 14b). The two illuminated edges of the wall are sep- 
arated by slightly less than a resolution element and, therefore, 


5.6. Constraints on the outer disk wall radius r 

Figure 15 shows the^ 2 plot for the outer disk wall radius r. The 
four epochs roughly agree with each other with a global best-fit 
value of r = 56 AU and well fitting solutions for the span of 
r = [52, 63] AU. This measurement is in excellent agreement 
with the value of r = 58 AU as derived solely on the basis of the 
SED by Espaillat et al. (2010). 
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Fig. 16. Constraints on the inclination i of the disk plane. The plot shows 
the excess of the best-fit p 2 for a given i with respect to the global 
minimum which is normalized by the threshold value Ap 2 . The 
color-coded curves represent each individual epoch, whereas the thick 
black curve is obtained by combining all four epochs. The bottom panel 
shows the well fitting range and the best-fit value for all five analyses. 


Fig. 17. Constraints on the Henyey-Greenstein parameter g. The plot 
shows the excess of the best-fitp 2 for a given g with respect to the global 
minimum p 2 in , which is normalized by the threshold value A;^ 2 . The 
color-coded curves represent each individual epoch, whereas the thick 
black curve is obtained by combining all four epochs. The bottom panel 
shows the well fitting range and the best-fit value for all five analyses. 


5.7. Constraints on the inclination i 

Figure 16 presents the x 2 plot f° r the inclination i of the disk 
plane. While there is some scatter among the lower-quality 
epochs, the combination of all four epochs yields a very well- 
behaved^ 2 profile with the best-fit solution at i = 50° and well 
fitting solutions spanning i = [44°, 54°]. This is in excellent 
agreement with the i = 51.5° ± 0.7° inferred by Pietu et al. 
(2007) from molecular line emission and confirmed by Andrews 
et al. (2011) but is notably higher than the i = 42° assumed by 
Espaillat et al. (2010). 

5.8. Constraints on the Henyey-Greenstein parameter g 

Figure 17 displays the p 2 plot for the Flenyey-Greenstein pa- 
rameter g, which describes the degree of anisotropic forward- 
scattering. There is a considerable amount of diversity among 
the epochs with local best-fits ranging from g = 0.60 to g = 0.86. 
Since g is a property of the size and shape distribution of the dust 
grains in the outer disk, it is not expected to exhibit any real as- 
trophysical variation. We therefore conclude that these variations 
represent measuring uncertainty and that the exact value of g is 
difficult to pinpoint with our analysis approach. Note that the 
Henyey-Greenstein formalism is a simplification of dust scat- 
tering behavior and may, thus, be an inaccurate description of 
the dust phase function at hand. We do not consider values of 
g > 0.9, since the Henyey-Greenstein formalism becomes very 
inaccurate in that regime. 

Nevertheless, combining the x 2 curves of the four epochs 
yields a clear best-fit solution at g = 0.67 and well fitting solu- 
tions for the range of g = [0.56,0.85]. All of these values are 
notably on the high end of the range expected for circumstellar 
disks, which explains the extremely high contrast observed be- 
tween the near and far sides of the disk in scattered light (see 
also Fig. 21). 


Because the system model is symmetric, one could argue that 
there is the possibility that the disk is oriented in the opposite 
way, where the far side being the bright side, and the grains have 
a negative assymmetry parameter g. The phase function of ag- 
gregate particles can be backward scattering over a limited range 
of scattering angles (see Min et al. 2010). In our case, the scatter- 
ing angles of the far side and the near side of the disk are approx- 
imately 134 and 34°, respectively (given that the opening angle 
of the disk at the location of the wall is approximately 6° and 
the inclination is 50°). To get the same contrast between these 
two scattering angles with a negative value of g, one would need 
g < -0.9, which is highly unphysical. 

5.9. Constraints on the center offset perpendicular to the line 
of nodes y 

Figure 18 shows the x 2 plot for the offset ij of the disk center 
from the star position perpendicular to the line of nodes. The 
four epochs agree well with each other. Combining thep 2 curves 
from all epochs yields a best-fit solution of y = -69 mas 
with well fitting solutions for a range of y = [-118, -44] mas. 
Accounting for the foreshortening due to the inclined disk, this 
implies a surprisingly large physical offset of r/ = 15 [10, 26] AU 
along the disk plane. Given the gap radius of r ~ 56 AU, 
this configuration corresponds to an eccentricity of e ~ 0.3. 
Figure 19 illustrates the difference between the eccentric best-fit 
solution for the K3 data and the restricted best-fit solution with 
y = 0 mas. 

Since this study is the first to measure the disk center offset 
along this dimension, its discovery is to be treated with caution 
until it can be confirmed independently. As Fig. 1 8 demonstrates, 
our combined^ 2 analysis rejects the default hypothesis of y = 0 
at the 7.5Ap 2 level; thus, there is little doubt that a strong y off- 
set is necessary for a best-fit solution within our model grid. 
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Fig. 18. Constraints on the disk center offset y perpendicular to the line 
of nodes. The plot shows the excess of the best-fit x 2 for a given y 
with respect to the global minimum X^ n - which is normalized by the 
threshold value Ax 2 - The color-coded curves represent each individual 
epoch, whereas the thick black curve is obtained by combining all four 
epochs. The bottom panel shows the well fitting range and the best-fit 
value for all five analyses. 




Fig. 19. Comparison of the best-fit solutions for y = -80 mas and y = 
0 mas for the K3 dataset, a) Image of the unconvolved best-fit model 
disk with y = -80 mas at logarithmic stretch, b) Same for y = 0 mas. 
c) S/N map of the residual image after subtracting the ADI-processed 
best-fit model image with y = -80 mas from the ADI-processed data 
at a linear stretch of ±5cr. d) Same for y = 0 mas. e) Difference of the 
images d)-c) at a linear stretch of ±2o\ The^ 2 of the y = 0 mas model 
exceeds that of the y = -80 mas model by 6.0 Ax 1 . 


However, we acknowledge the possibility that this offset repre- 
sents the model’s effort to match a feature of disk architecture 
beyond the scope of its parameter space. 

One degree of freedom that can conceivably be responsi- 
ble for an apparent shift of the disk image perpendicular to 
the line of nodes is the effective scale height of the outer disk, 
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represented in our model by the internal parameter Tsmau. As it 
increases, the upper and lower surfaces of the disk move apart in 
the projected image. Since only the upper surface is seen in our 
reflected-light images, a change in 4 , sma n may be perceived as a 
change in y in our model. 

The outer disk scale height is not treated as a free parame- 
ter of our model approach, since its default value v F sma ii = 0.5 
is constrained by the far-infrared SED (Mulders et al. 2010). 
However, we ran a small-scale parameter analysis to determine 
whether its introduction as a free model parameter would remove 
the need for a y offset. This analysis is described in Appendix C. 
To keep the nomenclature consistent with the other free model 
parameters, in particular, with the inner disk scale height s, we 
named the new free model parameter h := 2 v P sma n. This defines 
the default outer disk scale height as h = 1, which is analogous 
to the default inner disk scale height 5=1. 

In the small-scale parameter analysis, we explored the ex- 
treme case of h =2 (T' sma n = 1; that is, an inner disk scale 
height of twice the expected value). We found that the overall fit 
quality decreased significantly with respect to the default h = 1 
case, whereas the best-fit value of y did not change significantly 
from its h = 1 best-fit value. We therefore rejected the use of h 
as a free parameter and kept h = 1 fixed for our main analysis. A 
physical offset of the disk center from the star remains our best 
explanation for the observed disk images at this point. 

We note that our model implements an eccentric disk as a 
circular disk shifted away from the star. For low eccentricities 
e <k 1, this method approximates an elliptical gap well. At the 
observed eccentricity of e = 0.3, though, the approximation 
is quite inaccurate. We therefore acknowledge that our best-fit 
model likely cannot represent the disk’s true shape in all as- 
pects and that more elaborate disk models, including elliptical 
or even spiral-like disk architectures, may yield improved fits in 
the future. 

5.10. Covariances between r, i, g, w, y 

In Sect. 5.1, we opted for a brute-force parameter grid search 
in the five model dimensions r, i, g, w, y, since we expected 
that covariances between these parameters might render an inde- 
pendent optimization of each individual parameter impractical. 
Having calculated the x 2 values in this parameter grid for each 
epoch, we can now measure these covariances by mapping out 
the two-dimensional contours of the well fitting solution family 
for each pair of parameters. For a pair of independent parame- 
ters, the contours ideally take the shape of a laterally symmetric 
ellipse, whereas a positive [negative] covariance results in an el- 
lipse visibly skewed toward the rising [falling] diagonal. 

The results of this evaluation are presented in Fig. 20. Table 5 
lists the normalized correlation coefficients derived from the 
global well fitting solution family contours for each pair of pa- 
rameters. Note that the contours only comprise a limited number 
of grid points and are therefore only roughly accurate. 

We find the strongest covariance between (/, y). This positive 
covariance implies a negative covariance between (r,y), given 
that / := 56 AU/r. This behavior can be understood geometri- 
cally. The most prominent feature of the ADI images is the stark 
radial intensity gradient between the dark disk gap and the bright 
crescent of the forward-scattering near-side rim of the outer disk. 
As y is made more negative, the disk center is shifted further 
away from the star, bringing the disk’s bright near side closer to 
the star. To keep the radial position of the gap edge in its ob- 
served location, the disk radius r must increase to compensate. 
Decreasing the inclination i provides another way to widen the 
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Fig. 20. Visualization of the covariances between the model parameters r, i, g, w, y. The contours delineate the well fitting solution family 
(X 1 ^ X^i n + At 2 ) for each epoch and for the combination of all four epochs. The strongest covariances are seen between (/, y) and (y, g). 


Table 5. Correlation coefficients between pairs of model parameters in 
the brute-force grid-search study. 
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-0.55 


projected disk gap along the y axis and thus restore the position 
of the bright edge. Unlike a change in r, it changes the aspect 
ratio of the projected disk gap and thus degrades the fit to the 
curvature of the imaged bright crescent. 

The strong negative covariances between (/, g) and ( g , y) are 
more challenging to visualize. Decreasing y and / simultane- 
ously, as discussed in the previous paragraph, should result in an 
overall widening of the bright crescent in the ADI image, even 
though its shape and location is roughly kept constant. This may 
then require an increase in forward-scattering efficiency g to re- 
concentrate the flux in the central part of the crescent and thus 
restore the overall flux distribution. 

In principle, observing the ansae and the far side of the gap 
edge would break these degeneracies. However, those features 
are not discernible in our current datasets. Next-generation high- 
contrast imaging facilities like SPHERE ZIMPOL will likely be 


capable of such observations (Thalmann et al. 2008; de luan 
Ovelaret al. 2013). 

6. Discussion 

Overall, the disk architecture we derive from our NIR imag- 
ing data agrees well with the predictions from the SED (e.g., 
Espaillat et al. 2008) and millimeter and sub-millimeter interfer- 
ometry (e.g., Pietu et al. 2007; Andrews et al. 201 1). In the fol- 
lowing subsections, we focus on the disk aspects to which these 
previous methods of observations were insensitive, but which 
have become accessible through high-contrast imaging in re- 
flected light. 

6.1. Disk orientation (near vs. far side) 

While the position angle of the projected system axis is easily 
measured in millimeter interferometry (Pietu et al. 2007) and 
NIR imaging (Paper I), it is more challenging to determine which 
of side of the disk is inclined towards the viewer. Mulders et al. 
(2010) proposed two possible architectures for the LkCa 15 sys- 
tem, both of which are compatible with the SED constraints, 
but which predict opposite brightness asymmetries between the 
near and far sides of the disk. In one scenario, the bright cres- 
cent observed in NIR imaging represents the directly illuminated 
far side of the gap wall, whereas the near-side wall is obscured 
by the bulk of the optically thick disk. In the second scenario, 
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Fig. 21. Phase functions of dust grains discussed in this work. Small an- 
gles are forward-scattering; large angles are backward-scattering. The 
gray region corresponds to the best fit range of the Henyey-Greenstein 
asymmetry parameter g = [0.56,0.81]. The dotted line is the phase 
function of the dust grains used for the SED fit. The orange region in- 
dicates the observable range of angles from 34° to 1 34° . 


the observed crescent instead represents light forward-scattered 
off the near side of the outer disk rim, whose observable flux 
is greatly enhanced over the far side by anisotropic scattering. 
The far-side scenario is qualitatively supported by a disk mod- 
eling study by Jang-Condell & Turner (2013), whereas Pietu 
et al. (2007) and Grady et al. (2010) favor the near-side scenario 
based on asymmetries in the millimeter interferometry data and 
in space-based scattered-light images of the outer disk surface, 
respectively. Our first NIR imaging analysis in Paper I did not 
allow us to distinguish between these scenarios. 

The quantitative approach presented in this work, however, 
clearly favors the near-side scenario. We achieve visually and 
numerically convincing fits to the data using disk models with 
a high Henyey-Greenstein factor of g ~ 0.67, which results in 
an extreme enhancement of the near-side reflected light over the 
far-side flux. Conversely, while we were also able to create an 
opposite brightness asymmetry using the far-side scenario, most 
of the scattered light is constrained to a narrow, self-shadowing 
gap wall, which fails to reproduce the radially extended mor- 
phology of the flux in our images (Fig. 9). 

We therefore conclude that the northwestern side of the 
LkCa 15 disk is the near side, which agrees with Pietu et al. 
(2007) and Grady et al. (2010). 

6.2. Dust grain properties 

The derived value of the Henyey-Greenstein anisotropic scatter- 
ing parameter g ~ 0.67 [0.56,0.85] (cf. Fig. 17) can be used 
to estimate the size of the particles in the upper atmosphere 
of the disk, which are responsible for the scattering. Note that 
the inclination and opening angle of the LkCa 15 disk con- 
strain the observable phase function to the range of scattering 
angles roughly from 34° to 134° (Fig. 21). For larger particles, 
the anisotropy typically increases, but most of this effect is con- 
centrated at the smallest scattering angles (that is, the full phase 
function is no longer well represented with a Henyey-Greenstein 


function). At intermediate angles, the phase function usually flat- 
tens and becomes locally backward scattering for very large ag- 
gregates (negative values of g; see Min et al. 2010, and Min et al., 
in prep.). The rather large positive value of g we find for LkCa 15 
corresponds to particles that are roughly micron sized. This is in 
reasonable agreement with the grain sizes used in the upper at- 
mosphere of our model disk, and it implies that these grains are 
not agglomerated in large aggregates. 

6.3. Spatial asymmetry 

In Paper I, we postulated a physical offset x between the cen- 
ter of the LkCa 15 disk and the position of its host star along 
the line of nodes, based on the lateral asymmetry of the H - band 
ADI image. In our model, we included two free parameters (x, y) 
to allow an arbitrary displacement of the center of the circularly 
symmetric model disk from the star, where x is measured along 
the line of nodes (the major axis of the projected disk gap) and y 
perpendicular to it (along the minor axis of the projected disk 

gap)- 

Our analysis of the new A" s -band observations consistently 
yielded positive best-fit values of x (cf. Fig. 7), which tentatively 
confirm the existence of a lateral asymmetry in the appearance of 
the LkCa 15 disk. However, the absolute values of x show a con- 
siderable spread with the lower-quality epochs K2 and K4 sup- 
porting the large offset seen in Paper I, while the higher-quality 
epochs K1 and K3 advocate a much smaller offset consistent 
with zero. This may indicate that the real disk gap is roughly cen- 
tered along the x direction and that a feature further out on the 
disk surface is responsible for the observed asymmetry. For in- 
stance, a spiral density wave extending beyond the western ansa 
could shift the perceived center of gravity of the disk flux to- 
ward that side, biasing the analysis toward positive x values in 
the lower-quality datasets. 

In addition, our analysis suggests a previously unknown and 
suprisingly large offset in the y direction (cf. Fig. 18). Both 
the sign and the absolute value are consistent among all four 
l?s-band epochs. Correcting for the foreshortening along the in- 
clined disk plane, the best-fit offset of -69 mas translates to a 
physical displacement of ~15 AU, corresponding to ~0.3 times 
the gap radius. 

This is a dramatic departure from the roughly symmetric 
disk architecture that has been assumed in all previous studies 
of LkCa 15, and thus we treat this result with caution until it can 
be independently confirmed. However, we note that a sizeable 
disk asymmetry could well have gone undetected in previous ob- 
servations. Measurements of the system’s SED, as employed by 
Espaillat et al. (2008, 2010), are inherently insensitive to purely 
spatial information. While millimeter (Pietu et al. 2007) and sub- 
mm interferometry (Andrews et al. 2011) are capable of mea- 
suring disk information on large spatial scales at high S/Ns, the 
stellar photosphere is invisible at those wavelengths, rendering 
a direct comparison of the disk and star positions impossible. 
Andrews et al. (2011) do not attempt to fit any such offsets in 
their data but note that their disk center is <70 mas from the 
expected stellar position based on absolute pointing accuracy. 
This is marginally consistent with our best-fit offset. We also 
note that the near side of the disk appears brighter than the far 
side in Andrews et al. (2011), which is qualitatively expected if 
the near-side gap edge is closer to the star and, thus, has a higher 
equilibrium temperature. Since millimeter interferometry traces 
thermal emission rather than reflected light, anisotropic scatter- 
ing does account for this asymmetry. 
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Overall, our model treatment of the disk’s eccentricity - an 
azimuthally symmetric disk translated with respect to the star 
along the disk plane - is likely a simplification, whereas the real 
disk may have an elliptical gap, spiral arms, or even more irreg- 
ular features (cf. HD 142527; Canovas et al. 2013). However, 
the fact that in our analysis all four epochs of observation con- 
sistently suggests a significant offset is a strong indication of an 
underlying physical asymmetry in the system architecture, re- 
gardless of its detailed morphology. 

6.4. Wall shape 

To date and to our best knowledge, the shape of a transi- 
tional disk’s gap wall has only been detected in two other ob- 
jects: HD 100546 (Panic et al. 2014; Mulders et al. 2013b) and 
TW Hya (Ratzka et al. 2007, Menu et al., in prep.) with mid- 
infrared interferometry for both. In this work, we have inferred a 
wall-shape parameter of w ~ 0.30 for LkCa 15, which describes 
a gradual tapering of the outer disk surface density over a sig- 
nificant range of radii (a “round” or “fuzzy” wall) rather than a 
hard vertical cut-off ( w = 0). This value of w is comparable to 
the wall shapes described for the previously mentioned targets, 
indicating that similar processes are likely at work shaping these 
walls. It should be noted, however, that steeper wall shapes in 
other objects might have escaped detection as there may be less 
of a difference observationally with the commonly assumed ver- 
tical wall. On the other hand, the missing cavities in the SEEDS 
survey (Dong et al. 2012) might be interpreted as extremely ex- 
tended round walls. 

A tapered surface density profile for LkCa 15 was also 
proposed by Isella et al. (2009), as a way of reconciling the 
millimeter flux deficit with the existence of a NIR excess, with- 
out invoking an additional inner disk component. However, it 
requires an extremely round wall (w » 0.4) to achieve sufficient 
optical depth in the inner regions to create a NIR excess, which 
is not supported by our observations. Even though the disk wall 
extends over a wide radial range (cf. Fig. 5), the radiative trans- 
fer model requires an additional component to fit the NIR flux. 

Since these pre-transitional disks (transitional disks with a 
NIR excess, which implies the existence of an inner disk com- 
ponent at sub-AU radii) are thought to be shaped by planetary 
systems rather than photo-evaporation, it is tempting to explain 
the wall shape in terms of a planetary companion. Mulders et al. 
(2013b) have done so for HD 100546 using hydrodynamical 
models of planet-disk interactions and found that the extreme 
roundness ( w ~ 0.35) in this case is best explained by a brown 
dwarf companion. Taking their Fig. 7 at face value, the well fit- 
ting range of w = [0.19, 0.36] for LkCa 15 spans almost the en- 
tire plausible range of planet/star mass ratios 2 , making it difficult 
to establish constraints on planet mass. 

We note that planet mass also affects disk morphology in 
other ways, which we have not included in our model: 

Gcip depth. While our imaging data confirm an abrupt drop in 
surface brightness inside the LkCa 15 disk gap, they do not 
exclude the presence of a residual level of dust in that area. 
Knowing the degree of dust depletion in the gap (the “gap 
depth”) would impose additional constraints on the planet 
mass. However, the oversubtraction effects generated during 
the ADI data reduction by the stark gradients of the nearby 
bright gap render an absolute calibration of the flux level in 

2 One would need to take into account the difference in stellar mass be- 
tween LkCa 15 and HD 100546, which would scale down planet masses 
by a factor of ~2.5. 


the gap to be extremely difficult. Deep polarimetric imaging, 
on the other hand, may provide access to this information in 
the near future. 

Decoupling of gas and dust. In contrast to HD 100546, the 
(sub)micron-sized dust grains in the disk wall of LkCa 15 
show signs of dust settling (Espaillat et al. 2007; Mulders 
et al. 2010), indicating that they are dynamically decoupled 
from the gas and may be subject to radial drift, as in Pinilla 
et al. (2012b). This could alter the radial surface bright- 
ness/density profile of the disk and potentially provide a bet- 
ter diagnostic for planet mass. 

6.5. Constraints on planets and planet-induced dust grain 
differentiation 

As explored in Pinilla et al. (2012b) and de luan Ovelar et al. 
(2013), the presence of a planet in the disk generates a radial 
differentiation of dust grain sizes that is strongly related to the 
mass of the planet. The planet causes a pressure bump to form in 
the gas distribution of the disk, which traps large grains but al- 
lows for small grains to filter through towards smaller radii. As 
a result, the distribution of large grains exhibits a wider central 
gap than that of the small grains. Since NIR imaging and sub- 
mm interferometry mainly trace the micron- and mm-sized dust 
grains, respectively, this differentiation effect can be observed 
as a wavelength-dependent gap radius. Therefore, given the as- 
sumption that a single planet is responsible for the gap, the ratio 
of the gap radii measured at NIR and sub-mm wavelengths can 
be used to constrain the mass of the undetected planet (see Fig. 8 
in de Juan Ovelar et al. 2013). This scenario presents a plausi- 
ble explanation for the “missing cavities” sample of the SEEDS 
survey (Dong et al. 2012), where no cavities were found in polar- 
ized intensity //-band images of sources displaying large cavities 
in sub-mm studies. 

Our K s -band observations, the 870 pm observations pre- 
sented in Andrews et al. (201 1), and the SED fit carried out by 
Espaillat et al. (2010) all locate the wall at consistent radii of 
r ~ [50-58] AU. The object LkCa 15 is therefore different in 
this regard from the “missing cavities” sources. This lack of dust 
grain size differentiation strongly suggests that the radial differ- 
entiation of dust grain sizes, if present, is not significant. Figure 8 
of de Juan Ovelar et al. (2013) plots gap radius ratios between 
R-band and 880 pm as a function of planet mass. Preliminary 
calculations show that replacing R-band imaging with K s -band 
imaging does not significantly change this function. Taken at 
face value, gap radius ratios close to unity predict a mass of 
M < 1 Mj L | P for the planet responsible for the pressure distri- 
bution in the outer disc. In the case of a multi-planet system, this 
role is played by the outermost planet, as planets tend to influ- 
ence the disk locally (Goodman & Rafikov 2001). 

It may be worth it to note that the gap measured in the 
Andrews et al. (2011) study at 870 pm (~50 AU) is the smallest 
of all gaps measured at different wavelengths, which is interest- 
ing, since the substellar companion scenario assumed by both 
Pinilla et al. (2012b) and de Juan Ovelar et al. (2013) causes 
the opposite effect (that is, larger grains are further out). In their 
study, Andrews et al. (201 1) do mention that their uncertainties 
are underestimated, which could bring the value closer to the 
ones obtained in Espaillat et al. (2010) and this study; however, 
if confirmed, this effect would require further investigation be- 
fore any conclusions regarding the cause of this “reversed” dust 
grain size separation can be drawn. 

Simulations show that a single companion would need to be 
very massive to open a gap as large as the one observed in the 
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LkCa 15 disk ( M » 10 Mj up ; Kraus & Ireland 2012; Crida et al. 
2006). Such a companion would be at odds with the lack of ob- 
servable radial dust differentiation. However, a system compris- 
ing multiple low-mass planets could accommodate both of these 
constraints. Espaillat et al. (2008) note that the inner and outer 
boundaries of the gap in the LkCa 15 disk system are compara- 
ble to the smallest and largest orbital radii in our own planetary 
system, allowing for the tantalizing possibility that the LkCa 15 
system might be a young analog of our solar system. 

Kraus & Ireland (2012) propose a mass of ~6 M Jup for their 
planet candidate found in SAM observations. However, the au- 
thors note that this planet alone would be incapable of clear- 
ing the observed disk gap, which necessitates a second, lower- 
mass planet in a wider orbit closer to the disk gap. Since this 
outer planet would dominate the radial dust differentiation, our 
findings do not contradict the existence of the proposed massive 
planet. 

We cannot directly confirm or falsify the existence of the 
planet candidate reported by Kraus & Ireland (2012) from 
SAM observations in L' and K band, since its separation from 
the star (70-100 mas) lies within the inner working angle of our 
ADI observations. 

However, given the extreme anisotropy seen in our scattered- 
light images, the possibility should, perhaps, not be excluded 
that the disk itself may cause some of the structure seen in the 
SAM data, which has been inferred for the similar case of T Cha 
(Olofsson et al. 2013). 

6.6. Inner disk 

As first described in Mulders et al. (2010), the nature of the inner 
dust component of the LkCa 15 pre-transitional disk affects the 
observational appearance of the outer disk and can therefore be 
constrained by imaging observations of the outer disk. An op- 
tically thick inner disk would cast a shadow on the outer disk, 
whereas an optically thin dust shell would not. Espaillat et al. 
(201 1) report a time variability in the spectrum of the LkCa 15 
system in which the near- and mid-infrared fluxes exhibit neg- 
ative covariance (“see-saw effect”), which is taken as evidence 
for time-variable shadowing and, therefore, supports the opti- 
cally thick inner disk scenario. 

The effective scale height of the inner disk s is a free pa- 
rameter in our model; thus, our analysis is capable of detecting 
variable shadowing between epochs. Surprisingly,^ 2 minimiza- 
tion of this parameter (cf. Sect. 5.5) consistently yields s = 0 
for all epochs; that is, no shadowing is detected. This behavior 
can be explained with either the optically thin dust halo scenario 
(Mulders et al. 2010) or an optically thick inner disk tilted with 
respect to the outer disk, which causes its shadow to miss the 
outer disk wall at most azimuths. Both explanations are, how- 
ever, at odds with the “see-saw” behavior observed in the mid- 
infrared, which implies a photometrically significant amount of 
shadowing on the outer disk. Although s = 0 is a firm result of 
our analysis in all epochs, it is difficult to assess why the absence 
of a shadow increases the fit quality and which model parameters 
could alleviate the need for a fully illuminated outer disk. 

6. 7. Remaining uncertainties 

We caution that our forward-modeling analysis is limited by the 
extent of the disk model and its sampled parameter space. While 
we consider our best-fit models a convincing fit for our data, we 
acknowledge the possiblity that they may be approximations of 
a disk architecture beyond the scope of our model. 


Next-generation adaptive-optics facilities may offer the 
means to reveal the structure of the LkCa 15 disk gap with more 
conservative high-contrast imaging methods, resolving these re- 
maining uncertainties. Polarimetry is a particularly promising 
technique for such applications, as demonstrated by Quanz et al. 
(2013). 

7. Conclusion 

Following up on our H - band results from Paper I (epoch H 1 ), we 
have obtained four new epochs of A" s -band high-contrast direct- 
imaging observations of the LkCa 15 system to reveal new in- 
sights on the architecture of its pre-transitional disk. All four 
Ks-band datasets (epochs K1-K4) were taken in consistent ob- 
serving modes using Gemini NIRI with pupil tracking so as to 
enable direct comparison of the resulting ADI images. 

To compensate for the known flux loss and morphologi- 
cal impact of ADI, we have performed an extensive forward- 
modeling analysis. By simulating a parametric grid of scattered- 
light images of disk models, forward-modeling them through the 
ADI process, and comparing the results to the observed ADI im- 
ages with a x 2 metric, we have arrived at stringent quantitative 
constraints on the disk geometry. 

These results are independent from but consistent with the 
constraints previously inferred from the system’s SED (Espaillat 
et al. 2010) and from (sub-)millimeter interferometry (Andrews 
et al. 2011). However, they also include new findings that have 
only become accessible through NIR high-contrast imaging: 

- We find that the bright crescent seen in scattered-light ADI 
observations of LkCa 15 represents the strongly forward- 
scattering near side of the outer disk rather than the far side 
(cf. Sect. 5.2). 

- We tentatively confirm the existence of an asymmetry be- 
tween the ansae as postulated in Paper I. Assuming that the 
disk can be described as azimuthally symmetric but physi- 
cally displaced along the line of nodes (the major axis of the 
projected disk gap) by an offset x, we consistently find posi- 
tive best-fit values for x in all four A" s -band epochs, which is 
in line with our HI results (cf. Sect. 5.2). However, while the 
size of the best-fit x offset is consistent with the HI value in 
our lower-quality epochs K2 and K4, the higher-quality data 
of K1 and K3 indicate a much smaller offset consistent with 
zero. This may suggest that the observed asymmetry may be 
caused by a different mechanism than a displaced disk gap, 
such as, perhaps, by a spiral density wave beyond the west- 
ern ansa of the disk gap, which may be mistaken for the disk 
gap in lower-quality data. 

- Furthermore, our best-fit models exhibit a surprisingly large 
disk offset y that is perpendicular to the line of nodes (along 
the minor axis of the projected disk gap), bringing the 
brightly illuminated near-side rim of the outer disk closer to 
the star (cf. Sect. 5.9). Unlike x, the values of y are consistent 
with each other and significantly inconsistent with zero in all 
four Ks-band epochs, lending credibility to the result. While 
such a large displacement (~0.3 times the gap radius, depro- 
jected) constitutes a dramatic departure from the symmetric 
architecture commonly assumed for LkCa 15, it plausibly 
could have escaped detection in SED and sub-mm interfer- 
ometry. However, we note that a displaced circular gap is no 
longer a good approximation of an elliptical gap at an eccen- 
tricity of e = 0.3, so the real shape of the LkCa 15 disk may 
well lie outside of the model parameter space considered in 
this analysis. 
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- We find evidence for a “round” gap wall; that is, a gradual 
rather than abrupt gradient of surface density between the 
gap and the outer disk, which is consistent with the wall be- 
ing sculpted by a planetary companion (cf. Sect. 5.4). 

- We obtain best results for a completely absent shadow from 
the inner disk on the outer disk wall (cf. Sect. 5.5). This may 
indicate that the inner disk is inclined with respect to the 
outer disk, causing its shadow to miss the outer disk at most 
azimuths. 

The best-fit disk models are visualized in Fig. 10, and the nu- 
merical results of the analysis are summarized in Table 3. 

Finally, we acknowledge that the results of our forward- 
modeling analysis are bound by the scope of the model and 
its parameter space and that they are, therefore, simplistic ap- 
proximations to the real, complex reality of the LkCa 15 disk. 
Flowever, given the observational status quo, we consider them 
to represent the best estimate of the disk architecture to date. 

In the near future, imaging polarimetry with next-generation 
adaptive optics facilities, such as Subaru SCExAO (Guyon et al. 
201 1) or SPHERE ZIMPOL (Thalmann et al. 2008), may pro- 
vide a more detailed view of the LkCa 15 disk with fewer model 
dependencies. 


Acknowledgements. We thank F. Mem and S. Quanz for useful discussions, and 
the anonymous referee for helpful comments that improved the quality of the 
manuscript. C.Th. is supported by the European Commission under the Marie 
Curie IEF grant No. 329875. C.A.G. is supported by the US National Science 
Foundation under Award No. 1008440 and through the NASA Origins of Solar 
Systems program on NNG13PB64P. M.M. acknowledges funding from the EU 
FP7-2011 under Grant Agreement No. 284405. J.C. is supported by the US 
National Science Foundation under Award No. 1009203. This work is based 
on observations obtained at the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under a coopera- 
tive agreement with the NSF on behalf of the Gemini partnership: the National 
Science Foundation (United States), the National Research Council (Canada), 
CONICYT (Chile), the Australian Research Council (Australia), Ministerio da 
Ciencia, Tecnologia e Inova 5 ao (Brazil) and Ministerio de Ciencia, Tecnologia e 
Innovation Productiva (Argentina). The authors wish to recognize and acknowl- 
edge the very significant cultural role and reverence that the summit of Mauna 
Kea has always had within the indigenous Hawaiian community. We are most 
fortunate to have the opportunity to conduct observations from this mountain. 


References 

Amara, A., & Quanz, S. P. 2012, MNRAS, 427, 948 
Andrews, S. M„ & Williams, J. P. 2005, ApJ, 631, 1134 
Andrews, S. M., Rosenfeld, K. A., Wilner, D. J., & Bremer, M. 2011, ApJ, 742, 
L5 

Avenhaus, H., Quanz, S. P, Schmid, H. M., et al. 2014, ApJ, 781, 87 
Begemann, B., Dorschner, J., Henning, T., Mutschke, H., & Thamm, E. 1994, 
ApJ, 423, L71 

Bergin, E., Calvet, N., Sitko, M. L., et al. 2004, ApJ, 614, L133 
Boccaletti, A., Pantin, E., Lagrange, A.-M., et al. 2013, A&A, 560, A20 
Brandt, T. D., McElwain, M. W., Turner, E. L., et al. 2013, ApJ, 764, 183 
Buenzli, E., Thalmann, C., Vigan, A., et al. 2010, A&A, 524, LI 
Calvet, N., D’Alessio, P., Watson, D. M., et al. 2005, ApJ, 630, LI 85 
Canovas, H., Menard, F., Hales, A., et al. 2013, A&A, 556, A123 
Carson, J., Thalmann, C., Janson, M., et al. 2013, ApJ, 763, L32 
Crida, A., & Morbidelli, A. 2007, MNRAS, 377, 1324 
Crida, A., Morbidelli, A., & Masset, F. 2006, Icarus, 181, 587 
de Juan Ovelar, M., Min, M., Dominik, C., et al. 2013, A&A, 560, Al 1 1 
Dodson-Robinson, S. E., & Salyk, C. 2011, ApJ, 738, 131 
Dong, R., Rafikov, R., Zhu, Z., et al. 2012, ApJ, 750, 161 
Dorschner, J., Begemann, B., Henning, T., Jaeger, C., & Mutschke, H. 1995, 
A&A, 300, 503 


Duchene, G., McCabe, C., Ghez, A. M., & Macintosh, B. A. 2004, ApJ, 606, 
969 

Espaillat, C., Calvet, N., D’Alessio, P., et al. 2007, ApJ, 670, L135 
Espaillat, C., Calvet, N., Luhman, K. L., Muzerolle, J., & D’Alessio, P. 2008, 
ApJ, 682, LI 25 

Espaillat, C., D’Alessio, P., Hernandez, J., et al. 2010, ApJ, 717, 441 
Espaillat, C., Furlan, E., D’Alessio, P., et al. 2011, ApJ, 728, 49 
Fukagawa, M., Tamura, M., Itoh, Y., et al. 2010, PASJ, 62, 347 
Garufi, A., Quanz, S. P., Avenhaus, H., et al. 2013, A&A, 560, A 105 
Goodman, J., & Rafikov, R. R. 2001, ApJ, 552, 793 

Grady, C., Schneider, G., Sitko, M., & HST-GO-10177 Team 2010, in Proc. In 
the Spirit of Lyot 2010, ed. A. Bocaletti (Universite Paris Diderot) 

Guyon, O., Martinache, F., Clergeon, C., et al. 2011, Proc. SPIE, 8149, 08 
Henning, T., & Stognienko, R. 1996, A&A, 311, 291 
Henyey, L. G., & Greenstein, J. L. 1941, ApJ, 93, 70 
Hashimoto, J., Dong, R., Kudo, T., et al. 2012, ApJ, 758, L19 
Hodapp, K. W., Jensen, J. B., Irwin, E. M., et al. 2003, PASP, 1 15, 1388 
Huelamo, N., Lacour, S., Tuthill, P, et al. 2011, A&A, 528, L7 
Ingleby, L., Calvet, N., Herczeg, G., et al. 2013, ApJ, 767, 112 
Isella, A., Carpenter, J. M., & Sargent, A. I. 2009, ApJ, 701, 260 
Jang-Condell, H., & Turner, N. J. 2013, ApJ, 772, 34 
Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 
Kraus, A. L., & Ireland, M. J. 2012, ApJ, 745, 5 
Krijt, S., & Dominik, C. 2011, A&A, 531, A80 
Kuchner, M. J., & Holman, M. J. 2003, ApJ, 588, 1 1 10 
Lafreniere, D., Marois, C., Doyon, R., Nadeau, D., & Artigau, E. 2007, ApJ, 
660, 770 

Lagrange, A.-M., Bonnefoy, M., Chauvin, G., et al. 2010, Science, 329, 57 
Lubow, S. H., & D’Angelo, G. 2006, ApJ, 641, 526 

Luhman, K. L., Allen, P. R., Espaillat, C., Hartmann, L., & Calvet, N. 2010, 
ApJS, 186, 111 

Marois, C., Lafreniere, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 
641,556 

Marois, C., Macintosh, B., Barman, T., et al. 2008, Science, 322, 1348 
Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 
Mayama, S., Hashimoto, J., Muto, T., et al. 2012, ApJ, 760, L26 
Milli, J., Mouillet, D., Lagrange, A.-M., et al. 2012, A&A, 545, Alll 
Min, M., Hovenier, J. W., Waters, L. B. F. M., & de Koter, A. 2008, A&A, 489, 
135 

Min, M., Dullemond, C. R, Dominik, C., de Koter, A., & Hovenier, J. W. 2009, 
A&A, 497, 155 

Min, M., Kama, M., Dominik, C., & Waters, L. B. F. M. 2010, A&A, 509, L6 
Min, M., Dullemond, C. P., Kama, M., & Dominik, C. 2011, Icarus, 212, 416 
Min, M., Canovas, H., Mulders, G. D„ & Keller, C. U. 2012, A&A, 537, A75 
Mulders, G. D., Dominik, C., & Min, M. 2010, A&A, 512, Al 1 
Mulders, G. D., Min, M., Dominik, C., Debes, J. H., & Schneider, G. 2013a, 
A&A, 549, All 2 

Mulders, G. D., Paardekooper, S.-J., Panic, O., et al. 2013b, A&A, 557, A68 
Mutschke, H., Begemann, B., Dorschner, J., et al. 1998, A&A, 333, 188 
Olofsson, J., Benisty, M., Le Bouquin, J.-B., et al. 2013, A&A, 552, A4 
Panic, O., Ratzka, T., Mulders, G. D., et al. 2014, A&A, 562, A101 
Pinilla, R, Benisty, M., & Birnstiel, T. 2012, A&A, 545, A81 
Pietu, V., Dutrey, A., Guilloteau, S., Chapillon, E., & Pety, J. 2006, A&A, 460, 
L43 

Pietu, V., Dutrey, A., & Guilloteau, S. 2007, A&A, 467, 163 
Preibisch, T., Ossenkopf, V., Yorke, H. W., & Henning, T. 1993, A&A, 279, 577 
Quanz, S. R, Avenhaus, H., Buenzli, E., et al. 2013, ApJ, 766, L2 
Quillen, A. C. 2006, MNRAS, 372, L14 
Ratzka, T., Leinert, C., Henning, T., et al. 2007, A&A, 471, 173 
Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1 163 
Soummer, R., Pueyo, L., & Larkin, J. 2012, ApJ, 755, L28 
Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, M. F. 1989, AJ, 
97, 1451 

Suzuki, R., Kudo, T., Hashimoto, J., et al. 2010, Proc. SPIE, 7735, 30 
Thalmann, C., Schmid, H. M., Boccaletti, A., et al. 2008, Proc. SPIE, 7014, 3F 
Thalmann, C., Carson, J., Janson, M., et al. 2009, ApJ, 707, L123 
Thalmann, C., Grady, C. A., Goto, M., et al. 2010, ApJ, 718, L87 (Paper I) 
Thalmann, C., Janson, M., Buenzli, E., et al. 2011, ApJ, 743, L6 
Thalmann, C., Janson, M., Buenzli, E., et al. 2013, ApJ, 763, L29 
Van der Hulst, H. C. 1957, Light Scattering by Small Particles (New York: John 
Wiley & Sons, Inc.) 

Weaver, W. B., & Jones, G. 1992, ApJSS, 78, 239 

Zhu, Z., Nelson, R. R, Hartmann, L., Espaillat, C., & Calvet, N. 2011, ApJ, 729, 
47 


Pages 22 to 23 are available in the electronic edition of the journal at http://www.aanda.org 


A51, page 21 of 23 


A&A566, A51 (2014) 


Table A.l. Description of the preliminary best-fit (PBF) solutions to the 
four datasets. 


Parameter 

K1 

K2 

K3 

K4 

Scale factor / 

0.93 

1.00 

1.00 

0.97 

Gap radius r (AU) 

52 

56 

56 

54 

Inclination i (°) 

52 

50 

48 

48 

Henyey-Greenstein number g 

0.60 

0.80 

0.70 

0.60 

Settling parameter s 

0 

0 

0 

0 

Wall roundness w 

0.25 

0.35 

0.25 

0.30 

Orientation angle o 

147.7 

150.5 

150.7 

150.7 

Center offset x (mas) 

12 

73 

10 

45 

Center offset y (mas) 

-68 

-88 

-84 

-88 


Appendix A: Preliminary best-fit solutions 

The parameter sets listed in Table A. 1 represent the preliminary, 
coarse model fits to our four observing epochs of LkCa 15 ob- 
tained through an iterative search during the early phase of this 
work. They are used as a starting point for the studies presented 
in Sects. 5.2 and 5.5. Our final, optimized solutions and their 
confidence intervals are described in Sect. 5.3. 

Appendix B: Best-fit solutions under the restriction 
of y = 0 

Our analysis clearly favors a significant offset y of the disk center 
that is perpendicular to the line of nodes (cf. Sect. 5.9). However, 
this offset is yet unconfirmed, and could therefore represent an 
unknown bias in our modeling approach. For this reason, we 
here provide an overview of how the imposed restriction of y = 0 
would affect our results. 

Since some covariances are found between the parameters r, 
i, g, w and y (cf. Sect. 5.10), the^ 2 plots for those parameters 
undergo significant changes. Figures B.1-B.4 show the corre- 
sponding plots for the y = 0 parameter sub-space. 

In the unrestricted analysis, a range of y values are included 
in the well fitting family, each of which contributes a narrow 
valley to the x 2 landscape. For parameters that covary with y, 
the x 2 valleys change their positions with varying y; thus, the 
final x 2 valley calculated over all y is broadened. For this reason, 
imposing the restriction of y = 0 results in steeper^ 2 curves than 
in the unrestricted case. 

Furthermore, since the natural best-fit value of y is nega- 
tive, the restriction of y = 0 causes a negative [positive] shift 
in the best-fit value of parameters that exhibit positive [nega- 
tive] covariation with y. The gap radius r decreases from 56 AU 
to 50 AU, whereas the inclination i rises from 50° to 54°. 

Appendix C: Exploration of disk wall scale height h 
as an additional free parameter 

We ran a small-scale version of the brute-force parameter grid 
with outer disk scale heights of h = [1,2] to test whether the 
observed y offset could be an aliasing effect of an inaccurate 
assumption on h. Since h = 1 represents the hydrostatic equi- 
librium value, the test value of h = 2 is an extreme case that 
is intended to probe the possible codependence of h and y with 
high sensitivity. 

The analysis was run only for the best epoch, K3. We used 
the following parametric grid listed in Table C.l. Its irregularity 
occurs because it was assembled in two stages after it was found 
that the first stage did not fully encompass the best-fit solution. 


Scale factor f - 56 AU / r 
1.3 1.2 1.1 1.0 0.9 0.8 



Gap radius r (AU) 

Fig. B.l. Constraints on the outer disk wall radius r under the restriction 
of y = 0. The plot shows the excess of the best-fit x 2 f° r a given r 
with respect to the global minimum x ^ , which is normalized by the 
threshold value Ax 2 . The color-coded curves represent each individual 
epoch, whereas the thick black curve is obtained by combining all four 
epochs. The vertical dotted line marks the best-fit value of r from the 
unrestricted analysis. The bottom panel shows the well fitting range and 
the best-fit value for all five analyses. 



Inclination / (deg) 

Fig.B.2. Constraints on the inclination i of the disk plane under the re- 
striction of y = 0. The plot shows the excess of the best-fit x 2 for a 
given i with respect to the global minimum which is normalized 
by the threshold value Ax 2 . The color-coded curves represent each in- 
dividual epoch, whereas the thick black curve is obtained by combining 
all four epochs. The vertical dotted line marks the best-fit value of i 
from the unrestricted analysis. The bottom panel shows the well fitting 
range and the best-fit value for all five analyses. 
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Henyey-Greenstein g 

Fig.B.3. Constraints on the Henyey-Greenstein parameter g under the 
restriction of y = 0. The plot shows the excess of the best-fit y 2 for a 
given g with respect to the global minimum yT n , which is normalized 
by the threshold value Ay 2 . The color-coded curves represent each in- 
dividual epoch, whereas the thick black curve is obtained by combining 
all four epochs. The vertical dotted line marks the best-fit value of g 
from the unrestricted analysis. The bottom panel shows the well fitting 
range and the best-fit value for all five analyses. 



Wall shape parameter w 

Fig.B.4. Constraints on wall shape parameter w under the restriction 
of y = 0. The plot shows the excess of the best-fit y 2 for a given w 
with respect to the global minimum yT n , which is normalized by the 
threshold value Ay 2 . The color-coded curves represent each individual 
epoch, whereas the thick black curve is obtained by combining all four 
epochs. The vertical dotted line marks the best-fit value of w from the 
unrestricted analysis, which coincides with the values obtained under 
y = 0. The bottom panel shows the well fitting range and the best-fit 
value for all five analyses. 


Table C.l. Parameter grid for the h = 2 analysis. 


h = [1,21 

g = [0.3,0.5,0.7,0.8,0.91 

5 = [ 0 , 11 

w = [0.05,0.10,0.15,0.20,0.30] 

( = [40, 42, 44, 48, 52] deg 
/ = [0.933,0.967, 1.000, ..., 1.133] 
o = preliminary best-fit values 
x = preliminary best-fit values 

y = [-120,-110,-100,-90,-80,-60,-40,-20,0] mas 


The best-fit solution was achieved for g = 0.70, s = 0, w = 
0.05, i = 44°, / = 1.067, and y = -90 mas. The minimum y 2 
achieved was 601 , which is 2.8 Ay 2 above the y 1 . of the best-fit 
solution with h = 1 . While some of the best-fit parameters for 
h = 2 deviate significantly from their h = 1 counterparts (most 
notably the wall shape parameter w), we note that the best-fit 
value for y is still in excellent agreement with the results for 
h = 1 {y = -70 [-115,-45] mas). This implies that there is no 
significant covariance between y and /?. 

This behavior can be understood in conjunction with the ob- 
servation that the best-fit models do not include a shadow from 
the inner disk on the wall of the outer disk ( s = 0), and thus, 
the wall appears as a single bright crescent rather than two dis- 
tinct thin arcs. Increasing the outer disk scale height h makes 
the crescent wider (thus, perhaps, explaining the reduction in the 
wall shape parameter w, which is responsible for widening the 
crescent in the global best fit with h = 1 ) but does not affect its 
overall position. 

Thus, we conclude that the observed offset in y is a robust 
result of our analysis within the framework of our disk model 
and does not reflect an underlying error in the outer disk scale 
height h. We continue to use h = 1 as a fixed internal parameter 
for this work. 
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